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The theory of confinement and deconfinement is discussed as based on the properties 
of the QCD vacuum. The latter are described by field correlators of colour-electric and 
colour-magnetic fields in the vacuum, which can be calculated analytically and on the 
lattice. As a result one obtains a self-consistent theory of the confined region in the 
T) plane with realistic hadron properties. At the boundary of the confining region, 
the colour-electric confining correlator vanishes, and the remaining correlators describe 
strong nonperturbative dynamics in the deconfined region with (weakly) bound states. 
^ ' Resulting equation of state for /i = 0, p{T), are in good agreement with lattice 

' data. Phase transition occurs due to evaporation of a part of the colour-electric gluon 

, condensate, and the resulting critical temperatures Tc{tJ,) for different nf are in 

correspondence with available data. 
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1. Introduction 

The fundamental problems of confinement and deconfinement have recently become 
a hot topic because of a possible observation of deconfined phase — the Quark- 
Gluon Plasma (QGP) at RHIC The theory of the QGP was mainly centered 
on per turb ative ideas with an inclusion of some resummations and nonperturbative 
effects ^2l3] However, the experimental data^display a strong interacting QGP, with 
some properties more similar to a liquid, rather than to a weakly interacting plasma, 
with strong collective effects and large energy loss of fast constituents. These effects 
call for a nonperturbative treatment of QGP, and hence for a nonperturbative 
theory of confined and deconfined phases of QCD. 

For the confined phase the corresponding methods have been suggested in SI 
and formulated as the Field Correlator Method (FCM) — see ^ for a review. In 
the FCM at T = 0, the vacuum configurations support mostly two scalar quadratic 
in field correlators, D{x) and £'1(2:). The first correlator is purely nonperturbative 
and ensures confinement, while the second one contains both perturbative and 
nonperturbative parts. As a result one obtains linear confining potential from the 
D{x) plus perturbative and nonperturbative corrections coming from the Di{x). 
The spectra of all possible bound states (mesons, baryons, hybrids and glueballs) 
have been calculated in FCM, with current masses of quarks, the string tension a, 
and the strong coupling constant as(AQ) used as input, in good agreement with 
the experiment (see^^^for a review) and lattice data^. One can conclude therefore 
that the FCM describes well the confining region at T = /j, = 0. 

For T > several complications occur (see^^WH! for the details). First of all, 
colour-electric and colour-magnetic field correlators are no longer identical to each 
other and one has to do with five correlators instead of two: D^{x), Df{x), (x), 

(x), and Df^{x). Correspondingly, the colour-electric and colour- magnetic con- 
densates are now different, as well as the string tensions, which are simply related 
to the correlators. 



It is important to notice that only the electric tension a is responsible for con- 
finement, while all others ensure spin-orbit forces. Thus it was conjectured in ^ 
that deconfinement occurs if the correlator D^{x) vanishes at T = Tc{fJ,). Indeed, 
this conjecture was confirmed later on the lattice, while all other correlators were 
measured to stay intact Finally, it was argued in ^ (see also^ that, since the 
vacuum energy density is proportional to the gluonic condensate, one can derive 
the condition of vanishing of the D^{0) directly from the second thermodynamic 
law. Consequently, general features of the phase transitions and the very value of 
the critical temperature Tc{^) can be calculated in terms o f the difference of the 
gluonic condensates in the confining and deconfining phases EEl. 

We now come to the central point of the paper: namely what is the dynamics of 
the QGP and its Equation of State (EoS)? As will be shown below, the basic role 




(1) 
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in the EoS is played by two effects: the colour-electric (nonconfining) forces due 
to the correlator Df{x) and colour-magnetic forces due to the correlator D^{x). 
The former one is dominant for Tc < T < 1.5Tc and creates an effective mass 
(energy) of isolated quarks and gluons, thus leading to typical curves for the pres- 
sure p{T) similar to those observed in lattice calculations. Furthermore, the same 
Di{x) generates potential Vi{r) which is able to bind quarks and gluons ^^1^1^ and 
provides strong correlations in white Sq systems. Colour-magnetic forces which stem 
from the (x) also can bind quarks at large distances an d pr ovide a large ratio 
(potential energy)/ (kinetic energy) > 10 specific for liquids 1^^. Moreover, D^{x) 
(and thence cr^) grow with the temperature, and finally become the basic interac- 
tion in the diniensionally reduced 3-D QCD, with the transition temperature around 
1.5 2Tc. 

The paper is organised as follows. In Section[2]we introduce the FCM at T = 0. 
In particular, we give necessary essentials of the method, discuss colour-electric 
and colour-magnetic interactions, explain the QCD string approach to hadrons, 
and consider spin-dependent interactions. In Section [3] we generalise the FCM for 
nonzero temperatures: we introduce the winding measure of integration, define the 
Single Line Approximation (SLA), and derive the EoS of the QGP in the frame- 
work of the SLA. In Section [4] we discuss interactions between quarks and gluons 
above the Tc and investigate bound states of quarks and gluons due to these in- 
teractions. In Section [5] we generalise the results of Section [3] for nonzero chemical 
potentials /x > 0. In Section [6] we compare the results discussed in this review with 
other approaches found in the literature and discuss various solved and unsolved 
problems. 



2. FCM at T = 

2.1. Essentials of the method 

To describe the dynamics of quarks and gluons in both confining and deconfined 
phases of QCD we start from the gauge-covariant Green's function of a single quark 
(or gluon) in the field of other quarks and gluons plus vacuum fields and u se t he 
Fock-Feynman-Schwinger Representation (FFSR) (in Euclidean space-time) 

/•oo 

S{x,y\A) = {x\{m + D)-^\y) = {m-D) / ds{Dz),ye-'' ^a{x,y, s) 

Jo 

where K is the kinetic energy term, 

1 r . rdz^iry ^ 



(2) 



K = m^s + 



dr 



dT 



(3) 



with m being the pole mass of the quark. The parallel transporter along the tra- 
jectory 2p(t) of the quark propagating from point x to point y is given by 



^{x,y) = Pacxp 



dz^A^ 



(4) 
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(5) 



9 / <Jp-i^F^_,i,{z{T))dT 



and (T are the usual 2x2 Pauh matrices. The symbols Pa and Pf in ([4]) and ([5]) 
stand for the ordering operators for the matrices and Ffj_,y, respectively, along 
the quark path. 

Similarly, for one gluon Green's function, one writes 

G^,{x,y\A) = {x\{DlS^, - 2igF^,)-'\y) (7) 

and, proceeding in the same lines as for quarks, one arrives at the gluon Green's 
function in the FFSR: 

/•OO 

Gf,,{x,y\A)= / ds{Dz)^ye-''°^f,Ax,y,s), (8) 
Jo 

where 

^"^IJ^ [-^) '^^' ^^A^^y^s)^ ^x,y)PFexp(2g J dTF{z{T)) 

Here and below the tilde sign denotes quantities in the adjoint representation, for 
example, F^^l = iP-J-'^-. 

The single-quark(gluon) Green's functions ^ and ([5]) can be used now as build- 
ing blocks to write the Green's function of a hadron. For example, for the quark- 
antiquark meson one has: 

G,^{x, y\A) = Tr{S,{x, y\A)TSl{x, y\A)T^), (10) 

where F is a Dirac 4x4 matrix which provides the correct quantum numbers of 
the meson (F^l, 7^, 75, 7^75, . . .) and the trace is taken in both Dirac and colour 
indices. The next important step is building the physical Green's function of the 
meson, 

Gqqix, y) = (Gqg(x, ^1^4))^, (11) 

where the averaging over gluonic fields is done with the usual Euclidean weight con- 
taining all gauge-fixing and ghost terms. The exact form of this weight is inessential 
for what follows. 

The actual average one needs to evaluate in order to proceed reads: 

(C)) = ($(C)$,, {C, Si)$,, {G, S2))a, (12) 

where the closed contour G runs along the trajectories of the quark, zi^(ti), and 
that of the antiquark, Z2i^{t2). Since the orderings Pa and Pf in ^^i, and ^^2, are 
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universal, then Wcria2{C) is just a Wegner- Wilson loop with the insertion of the 
operators {criF) and {<T2F) at proper places along the contour C. 
To proceed it is convenient to rewrite (jl2p in the form 



{W,a'{C)) = (^expig j^dTi^,{z)F^,{z)J , (13) 

where the non- Abelian Stokes theorem was used and the integral is taken over 
the surface S bounded by the contour C. Notice that the differential 

d7r^^(z) ^ ds^,^{z) - ia'i^JdTi + ia^^}dT2 (14) 

contains not only the surface element ds^,y but it also incorporates the spin vari- 
ables. This is the most economical way to include into consideration spin-dependent 
interactions between quarks (ti and T2 are the proper-time variables for the quark 
and the antiquark, respectively). Furthermore, if a cluster expansion theorem is 
applied to the right-hand side of ([13]), then the result reads: 

^{W.^.AC)) = expf; j d^{l) ... 1 d^(n)((F(l) . . . F{n))), (15) 

where d7r(fc) = dn f_,^i,^{z) andi^(fc) = i^^^,,^^ (ufe, xq) = $(xo, Ufc)F^fc,.jMfc)$(ufc, xq). 
Double brackets ((. . .)) denote irreducible (connected) averages, and the reference 
point is arbitrary. 

Equation (|15p is exact and therefore its right-hand side does not depend on any 
particular choice of the surface S. At this step one can make the first approximation 
by keeping only the lowest (quadratic, or Gaussian) field correlator {{FF)\ while 
the surface is chosen to be the minimal area surface. As it was argued in^iSl, using 
comparison with lattice data, this approximation (sometimes called the Gaussian 
approximation) has the accuracy of a few per cent. 

As was mentioned before, the factors (m — D) in the Green's function pT|l need 
a spec ial treatment in the process of averaging over the gluonic fields — it is shown 
in that one can use a simple replacement, 

m-D^m-ip, Pt^ ^ \ {~^^ ■ '^^^^ 

Therefore, the first step is fulfilled, namely the derivation of the physical Green's 
function of a quark-antiquark meson (one can proceed along the same lines for 
baryons and hadrons with an excited glue) in terms of the vacuum correlator 
{{FF)). In the next chapter the structure of this Gaussian correlator is studied, 
in particular, its separation into colour-electric and colour-magnetic parts. 
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2.2. Colour- electric and colour-magnetic correlators 



The Gaussian correlator of background gluonic fields can be parameterized via two 
scalar functions, D and Di 31 

9' 



{{FF)) EE D,.xp{z - z') 



|p((TrF,.(z)$(z,z')i^Ap(^')$n^,^')» 



{Sf_t\S^p - Sf,pS^x)D{u) + i 



d 



A 



(17) 
Diiu), 



^ — (uxSiyp - UpSxu) + 
OUfi \A ^ p y 

where u = z — z' . Notice that, in order to proceed from (fT5|) to (flTll one needs 
to show that, for the correlators in taken at close points, \z — z'\ < tJ^, the 
parallel transporters $ passing through the point xo can be reduced to the straight 
line between the points z and z' . To this end, notice that, for a generic configuration 
with \z — z'\ <^ \z — xq\,\z' — xa\ ^ R (i? is the radius of the Wilson loop, of order 
of the hadron size), one can write $(z — xq)(^{xq — z') w $(z — z') + 0{Tg /R^). 
For future convenience let us also write the correlator (|17p in components: 



|-((Tri?,(z)ci>i?,(z')<i>^))-4 

9 

i{TTH,{z)<^H,{z')<^^)) = 5,, 



^{(YtH,{z)^E,{z')<^^)) 



D^{u) 



,dDf 



u 



UiUj 



UiUj- 



dv? ' 



(18) 



EH 



For the sake of brevity we omit, in this chapter, the spin-dependent terms in the 
Wilson loop (|15p and consider only nonperturbative contributions to the Gaussian 
correlator {{FF)) responsible for confinement and given by the correlator D. This 
correlator enters the Wilson loop multiplied by the surface elements dsp^ds\p 
and, in what follows, we distinguish between the colour-electric and colour-magnetic 
contributions in psp . The former are accompanied by the structure dsoidsQi and 
enters (ITSI) multiplied by electric correlator , whereas, for the latter, these are 
dsjkdsjk and , respectively. The electric and magnetic string tensions are defined 
then according to ([1]). 

We now proceed from the proper-time variables to the laboratory times by 
performing the following change of variables: 

dziQ J dz2o 



dri 



dT2 



(19) 



2^1 2^2 

and synchronise the quarks in the laboratory frame, putting ziq — Z20 — t. The 
new variables Ml. 2 which appeared in (jl9p are known as the auxiliary (or einbein) 
fields The physical meaning and the role played by the einbeins will be 



'^The gluonic correlation length Tg defines the distance at which the background gluonic fields are 
correlated or, more specifically, the correlator D{u) defined in I I17II decreases in all directions of 
the Euclidean space, and the length Tg governs t his dec rease. The correlation length Tg extracted 
from the lattice data is quite small, Tg < 0.1 fml ^S | 19 | 
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discussed in detail below. Now we can write the surface element through the string 
profile function w^{t,f3) as 

ds^,{z) = e'^^'daw^it, P)dbW,it, (3)dtd(3, {a, b} = {t, /?}, (20) 

where ^ t ^ imax, 0^/3 ^ 1, and we adopt the straight-line ansatz for the 
minimal string profile, so that the latter is defined by the trajectories of the quarks 

US 

Wf,it, 13) = (izi^it) + (1 - P)z2^,{t). (21) 
For further convenience let us introduce two vectors: 

r^Zi-Z2, p = [(^1 - ^2) X (/3i;i + (1 - /3)i2)] = rw, (22) 

where u) is the angular rotation vector. This allows one to write the differentials in 
a compact form, 

dsoi{z)dsQi{z') = r{t)r{t')dtdt' dfidp' , ds.jk{z)dsjk{z') = 2p{t)p{t')dtdt' d(3d(3' . 

(23) 

Presenting the averaged Wilson loop (fT5|) (without spins) as 
^{TrW{C))=e-', J ^ + , 
one can write for the electric and magnetic contributions separately: 

J^=/" dtdt' [ dp d(3' r{t)r{t')D^[z- z'), (24) 
Jo Jo 

^ f '^^'^ dt dt' f dp df3' p{t,l3)p{t\(3')D"{z- z'). (25) 



The correlation functions D^'^ decrease in all directions of the Euclidean space- 
time with the correlation length Tg which is measure d on th e lattice to be rather 
smaU, Tg « 0.2-=- 0.3 fmlSlor even smaller, Tg < 0.1 fmHSHH. Therefore, only close 
points z and z' are correlated, so that one can neglect the difference between r{t) 
and r{t'), pit, (3) and p{t',P') in ^ and also write: 

{z-z'f^{z{t,p)-z{t',(3')f=g'^'Ub, ia^t~t', 6 = (26) 

The induced metric tensor is g"'' = g°-5°''', g^g^ — det g — + p'^ = r^(l -I- lu^). 
Now, after an appropriate change of variables and introducing the string tensions, 
according to one readily finds: 

= dt f dp jl-^ = aET dt C dp^=L=, (27) 

Jo Jo + p2 Jo Jo VI + uj^ 

J" ^an I dt I dp ^4__ = anr / dt / dp (28) 



y/r^ + f^ Jo Jo vl 
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For as = cth — c the sum of and repr oduces the well-known action of 
the Nambu-Goto string (see, for example, ^201-1221)^ 



J^J^ + = ar / dt df3y/l+uj^. (29) 
Jo Jo 

Notice that, for oj 1, Je = (Jet x ty^ax + O(w^) and Jh = 0{uj'^), that is con- 
finement is of a purely colour-electric nature while the colour-magnetic contribution 
is entirely due to the rotation of the system. 



2.3. QCD string approach 



With the form of the colour-electric and colour-magnetic spin-independent inter- 
actions (|27p and ([28]) in hand we are in a position to build a quantum-mechanical 
model of hadron consisting of quarks (gluons) connected by straight-line Nambu- 
Goto strings. Below we consider, as a paradigmatic quark-antiquark meson. 
We start from the kinetic energies of the quarks, 

2 



TO,- Si + 



1 



dTi 



dTi 



rni_ iM 
2 



dt '^.^^^^ ^ i = \X 

(30) 

where (|19|) was used, and path integrals in /ii appear through the substitution 
dsiDzio —^ DfjLi. One can see therefore that the einbeins should be treated on equal 
footings with other coordinates. However, since time derivatives of the einbeins, /ti, 
do not enter ([50]) , they can be eliminated with the help of the well-known integral: 



ZJ/i exp 



exp 



ah dt 



(31) 



Then, with the help of ([2]), ((27|) . ([28]), and ([30]) . one can extract the Lagrangian 
of the quark-antiquark system in the form (notice that hereinafter in this chapter 
we work in Minkowski space-time); 



L 



TOiW 1 - 2^ - m2 



1 — z\ — Get 



[n X (/3ii + (1 - /3)Z2)]^ 



df3- 



(32) 



y/l-[nx iPzi + (1 - P)Z2W 
In fac t, the einbein form of the kinetic energies is much more convenient 



[201- 1221231 



since it does not contain square roots (unbearable for path integrals) 



and one deals with formally nonrelativistic kinetic terms, while the entire set of 
relativistic corrections is summed by taking extrema in the einbeins. Furthermore, 
extra (continuous) einbeins, viP) and r]{P), can be introduced in order to simplify 
the string terms in ([32]) . through the substitutions: 




d(3- 



^4! 

B 



d(3 {Brf + 2Ari) . (33) 
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Then, introducing the centre-of-mass position and the relative coordinate as 

/^i + fn I3vdd , ^ 

H = Cia^i + (l-Ci)a^2, r = Xi-X2, Ci^— ,i , C2 = l^Ci, (34) 

Ml + M2 + Jo i^dp 

one can rewrite the Lagrangian in the form (the centre-of-mass motion is omitted 
for simpUcity): 

^ - g [g + I] - /' + f + + jM^n)^ + l^i- X nf. ,35) 

where ai = an +'I]'^{(Jh — o'e), 0^2 — 2?7(ct£; — an) and we have defined the reduced 
masses for the angular and radial motion separately: 

M= \ , 36) 

^J■l + M2 



M = A^i(l-Ci)'+M2Cr+ / iP-Cir'^df]. (37) 

Jo 

The original form of the Lagrangian is readily restored once extrema in all four ein- 
beins, {/xi, /i2, i^(/?), r?(/3)}, are taken. Generally speaking, the einbein fields appear 
in the Lagrangian and, as was mentioned before, even in absence of the corre- 
sponding velocities, they can be considered as extra degrees of freedom introduced 
to the system. The einbeins can be touched upon when proceeding from the La- 
grangian of the system to its Hamiltonian and thus they mix with the ordinary 
particles coordinates and momenta. Besides, in order to preserve the number of 
physical degrees of freedom, constrains are to be imposed on the system and then 
the formalism of constrained systems quantisation is operative (see, for exam- 
ple, for the open straight -line QCD string quantisation using this formalism) . 
A nontrivial algebra of constraints and the process of disentangling the physical 
degrees of freedom and non-physical ones make the problem very complicated. In 
the meantime, a simpler approach to einbeins exists which amounts to considering 
all (or some) of them as variational parameters and thus to taking extrema in the 
einbeins either in the Hamiltonian or in its spectrum I 23 | 25 | ggjjjg an approximate 
approach this technique appears accurate enough (see, for example, providing 
a simple but powerful and intuitive method of investigation. The physical meaning 
of the variables /i^ (i = 1, 2) is the average kinetic energy o f the i- th particle in the 
given state, namely, iii = (i/p^ -|- m^) (see the discussion inl2I126l)_ xhe continuous 
einbein variable ^{(3) has the meaning of the QCD string energy density 

Following the standard procedure, we build now the canonical momentum as 

dL 

p = — ^ ^iinr)n + ji{r - n{nr)) , (38) 
or 

with its radial component and transverse component being 



(np) = fi{nr), [n x p] = /i[n x r], 



(39) 
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respectively. Thus we arrive at the spin-independent part of the Hamihonian l^H: 



2^j 



2 



d(3 



2v 



PI 
2n 



2/ir2 



(40) 



Notice that the kinetic part of the Hamihonian (|40p has a very clear structure: 
the radial motion of the quarks happens with the effective mass /z, whereas for the 
orbital motion the mass (1 is somewhat different, containing the contribution of the 
inertia of the string. For as = cth = o", the field r/ drops from the Hamiltonian and 

91 

the standard expression for the string with quarks at the ends readily comes out 
from (I40l). 



2.4. Spin- dependent interactions in hadrons 

In heavy quarko nia (up t o the order 1/m^) the standard Eichten-Feinberg decom- 
position is valid |27 | 28 | 29 j. 



4m2 



r dr 



(3(<Tin)(<T2n) - (Tio-a) 
12mim2 



2dVi 
r dr 



(T2)L 1 dV2 

2miTO2 r dr 



CT\<Jl 



(41) 



where each potential V^(r) (n =0-4) contains both perturbative (P) and nonper- 
turbative (NP) contributions: Vnif) = VJ^(?') + V^^{r'). The static interquark 
potential Vb (r) , together with the potentials V\ (r) and V2 (?') , satisfies the Gromes 
relation 



<{r) + V[(r)-V!,(r)^Q. 



(42) 



Notice that this relation refers both to the perturbative and nonperturbative parts 
of the potentials Vnij") (jt- = 0, 1, 2). 

We now return to the Green's function of the quark-antiquark system ([2]) and 
restore spin-dependent terms in order to derive the generalisation of (|4ip . 

Let us quote here without derivation the full set of spin-dependent potentials, 
both perturbative and nonperturbative, obtained in the framework of FCM and 
with the string rotation taken into account (an interested reader can find the details 
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of the derivation in'^H): 



(T-\ L (T oL 

— = 1 1— 



r dr 

X 

2dVi 
r dr 

X 



1 

2^ 



oo 

dv I dX 

Jo 



(i-Ci)^ + (i 
filfi 

dv 



D + Di 
(T2L 



C2) 







M2M 
d\D \ 1 



(i-Ci)^ + (i 



MlM 



C2) 



<j->L 



1 -1- aj)]LldV2 
r dr 



2mim2 r dr r 
(3(crin)(cr2n) - (Ti(T2 



\l dv 





XdX 



D 



M2M 
£>! + A 



2^^ 



aA2 



(q-i +<T2)X 



(TiCr2 



12miTO2 



12TOim2 
1^4 (r) ^ 



t/aW ^ -2r 



/>oo 

6 / diy 
Jo 



D{r, v) + 



dvDi{r, v) 


2 2_5_ 

3'' 



(3((Tin)(cr2n) - (TiCTa) 



12/11/^2 
(TiCr2 



12/ii^2 ■ 



(43) 

where the masses m.i are replaced by Hi andu, which makes this result applicable 
also to light quarks. Notice that this r esultM is not due to the 1/m expansion, but 
is obtained with the only approximation made being the Gaussian approximation 
for field correlators (corrections may come from triple an d qu artic correlato rs). 
Accuracy of this approximation was checked both at T = and at T > 
to be of the order of one percent. 

With this explicit form of the potentials and taking the limit of heavy quarks 
one can check the Gromes relation (|42p . which now reads: 



/•oo /*r 

V^{r) + V;{r) - V^{r) = 2 dv dX[D^{X, v) {X, v)] 

Jo Jo 



+r / dv[Df{r,v)-D^(r,v)]. 



(44) 



At T = 0, Z?^ = and Df = , so that the Gromes relation (|42]) is satisfied. 
In the next Section we shall discuss the FCM at nonzero temperatures, in particular 
at T > Tc 7 where the latter equalities between colour-electric and colour-magnetic 
correlators do not hold, and the Gromes relation is therefore violated. 



2, 2009 16:39 WSPC/INSTRUCTION FILE QGP 



Deconfinement and quark-gluon plasma 13 

3. FCM at T > 

3.1. The winding measure of integration 

Now we turn to the case T > and use Matsubara technique for the path integrals 
in the FFSR, first introduced inH^, 

N 



{Dzr^y = hm n 



1—1 ^ ^ 



ri=0,±l,... ^ ' 



(45) 



where C(fc) = z{k) — z{k — 1) and e = t/N , and (3 — 1/T. It is easy to check that, 
with such a measure, for example, the free massless propagator takes its standard 
form: 



dt exp 



TV 

E 



Ae 



n dCjm) d^p 



(47re)2 ^ (27r)4 
exp \ip Cim) - (x-y)- n/3(5^4^ 



(46) 



f d^\ 
j exp \_~p'^t — ip(a; — z) — ipin(i\ dtj^ 



d^p 



E 



Td'^p exp[— ipi(a; — y)i — i2iTkT(xi — yn)] 



(27r)3 



pj + (27rfcr)2 



fc=o,±i,. 

where the Poisson summation formula was used to obtain: 



exp(ip4n/3) = ^ 27r(5(p4/3 — 27rfc). 

ra=0,±l,... fc=0,±l,... 



(47) 



For quarks one is to take into account the antisymmetric nature of the fermionic 
fields, which yields instead of (|45)) 



n 



E 



n=0,±l,... 



(27r)4 



exp 



N 



ip j C(™) ~ (2; — y) — 



>4 



(48) 



Therefore, in order to proceed from T = to T > we simply replace {Dz)xy 

w 

by {Dz)'^y or {Dz)^^ in all expressions. 

We now come to the point where we need to separate the quantum giuon field 
ttfx from the vacuum background field B^, so that = + a^, both satisfying 
periodic boundary conditions, 



Bfj,{zi, z) = B^{zi + n/3, 2;), a^(z4, z) = a^{zi + n/3, 2), 



(49) 
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where n is an integer. The partition function can be written as 
Z{V,T) = Z{B)=N j D{<j))e^^ ("/""'^ j d?xUot{x,T)^ (50) 

where is a normaUzation constant and denotes the entire set of fields a^t, 5*, 
. The expUcit form of the Lagrangian Ltot is given in . Furthermore, in (|50p 

{■ ■ ■) B stands for the averaging over the (nonperturbative) background fields _B^. 

For our purposes one needs not know the details of this averaging. 

Integration over the ghost and gluon degrees of freedom in (|50|) yields the same 

answer as for the case of T = 0, but now all the fields are subject to the periodic 

boundary conditions (j49|) . Disregarding, for the sake of simplicity, quark and source 

terms in Ltot, one obtains 

Z{B) = N'{detWiB));^\/^[deti-D^{B)D4B + a))]a=s/sj 
^ ^ 1 . / SV) ( \ 



where the valence-gluon Green's function G = is G^jy — ■ 1 + 2igFfj_^)^^ . 
We can consider strong background fields, so that gi?^ is large (as compared to 
Aq(^£j), while as = /{Att) in this strong background is small at all distances. 
Moreover, it was shown that Us is frozen at large distances In this case (jSTj) is 
a perturbative sum in powers of g", arising from the expansion in (t/a^)". 

In what follows we shall discuss the Feynman graphs for the thermodynamic 
potential F{T,iJ.), related to Z{B) via 

F{T,fi) = -T\n{Z{B))B. (52) 



3.2. Nonperturbative dynamics of quarks and gluons in the 
single-line approximation 

In this chapter we come to explicit calculations of various properties of the quark- 
gluon plasma. It is clear that ([5^ contains all possible interactions, perturbative 
and nonperturbative, between quarks and gluons and, in particular, creation and 
dissociation of bound states. It is impossible to take into account all possible sub- 
systems and interaction between them, so it is imperative to choose the strategy 
of approximations for the quark-gluon plasma as a deconfined state of quarks and 
gluons. 

We assume that the whole system of Ng gluons, Nq quarks, and Nq antiquarks 
remains gauge-invariant at T > Tc, as it was at T < Tc. However, in case of decon- 
finement and neglecting (in the first approximation) all perturbative and nonper- 
turbative interactions, any white system possesses the same energy, depending only 
on the number and type of constituents. In this case the partition function factorises 
into a product of one-gluon and one-quark (antiquark) contributions, and one can 
calculate the corresponding thermodynamic potential. This first step is called the 
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Single Line Approximation (SLA) and below we shall calculate the results for the 
gas of quarks and eluons in the strong vacuum fields using our path-integral for- 
mahsm. We follow ^121351-137] 

In this chapter we consider gluons and quarks in SLA. Notice that, although 
the interaction between quarks and gluons is switched off in this approximation, 
there exists a strong interaction of quarks and gluons with the nonperturbative 
vacuum fields. We therefore keep only the lowest order in {ga^) but take into account 
interaction with the nonperturbative background vacuum fields to all orders, using 
the field correlator technique. As a result one can write 

Zo = e-^o(^)/^ = N'{eM-Fo{B)lT))B. (53) 

Then single-gluon and single-quark contributions to the free energy can be written 
in the form: 

^F^\B) = ilndctG"^ -lndct(-i:i2(B)) = 

= Tr {-1 [ e(.)^e--- -f e(.)^e-^(-)} (54) 

and 

^F^{B') ^\\i.Aet{ml~b\B')) = -ixr^" e(^)ye— (55) 

respectively, where i?^ = 3^ — g)S^4, fi is the quark chemical potential, and 
is the current quark mass (the pole mass when next terms in as are considered). 
The operator D^{B') can be presented as: 

D^iB') = (D^j^f = Dl{B') - gF^,a^, = - gaF. (56) 



Then, using the FFSR, one can write 
{F^'iB))B = -tJ ^^{s)d''xiDz):,e-'< 



(57) 

where taking trace implies summation over the Lorentz and colour indices and 
is defined in ([9]). 

Similarly, one has for the quarks: 

rri />00 7 

{F^iB'))B^-- -C(s)d4a.(D^le-^-^™'Tr(<i>,(x,x,s))s, (58) 
^ Jo * 

where Tr implies summation over the spin and colour indices and the operator 
is given in (O. 

For future convenience let us define the proper-time gluon and quark kernels, 
Gis) = JiDzre-''TT{Ws{C,s)), S{s) = J JUzf" e'^'TriW^iC s)) (59) 
where 



(N^ - 1) 



-(bp - $ 
2 



W.^—<i>.. (60) 



2, 2009 16:39 WSPC/INSTRUCTION FILE QGP 



16 A.V. Nefediev, Yu.A. Simonov, M.A. Trusov 

Then, using the relation between the pressure and the free energy, PV3 — 
— {Fq{B))b, one can find (the T-independent term with n = is subtracted, V3 
is set equal to unity) 

p,,=(7v,^-i) r^5:G(")(.) (61) 



for gluons and 

= 2Nc I —e-<' V(-l)"+M'5'"Hs) + S'^'"\s% (62) 



=1 



for quarks, respectively. 

Therefore, in or der to calculate G("Hs) and S'^''\s), one needs to evaluate the 
Wilson loop operators (W^s,o-(C'„)). First, we neglect the spin-dependent terms in 
these Wilson loops since spin-dependent interactions at T = were considered 
before, while their generalisation for T > is straightforward. Thus we deal with 
the operator {W{C)). Let us look more carefully into the topology of this Wilson 
loop. It is important to bear in mind that the basis states |fc) with the quark-gluon 
numbers Ng''^ , iV^'^'' , iV^'^-' entering the partition function. 



n\e-"/^\n), (63) 



are gauge-invariant, and one can use a generic decomposition of the type: \k) — 
\{99){99){999)i.1^) ■ ■ ■)^ where particles in parentheses form white combinations. 
Then, neglecting interactions between white subsystems, one finds that the contri- 
bution of such white subsystems, for example, {gg) and {qq), to be 



Z(i2) = j dridr2(TrI^(C(i\c(2))), dV, = ds,Dz,e- 



Ki 



where W{C^n^C^n^) is a closed Wilson loop made of paths of the two gluons, of the 
qq pair, and, finally, of the parallel transporters (Schwinger lines) in the initial and 
final states — the latter are necessary to make the states gauge-invariant. 

To proc eed we apply the non-Abelian Stokes theorem and the Gaussian approxi- 
mation (see^I^^l^for the discussion and relevant references) to arrive at (fT5)) . and it 
is important now to specify the surface Sn with the surface elements ds^^ (u) in the 
integral on the r.h.s. of (fT5|). For this analysis let us to consider colour-electric and 
colour-magnetic contributions separately. For the field correlator we use formulae 

m- 

Consider first the colour-electric piece with v — a ~ A (note that, by definition, 
fjL < V and X < p): 

-^n =7; Diiki{u - v)daii{u)dakA{v). (64) 



2 . 

The surface Sn is inside the winding Wilson loop for the gauge-invariant qq 
system. Since the correlator D^{u) vanishes at T > Tc, while the contribution 
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of the correlator Df leads, at large interquark distances, to a sum of terms for 
the quarks and antiquarks separately (we shall study the interaction of quarks in 
chapter m below). Then for a single quark one obtains in this way 

1 rnP rnf3 poo 



2 _ 
1 



n/3 



^np dv (^1 - ^ ^d^DfiVi^'+e)- (65) 

Note that, for n = 1, one recovers the expression for the Polyakov loop obtained 
in I35i37 ^ namely 

ifund = exp (^-l,yi(oo)^ , 2^1"i(oo) = Jf . (66) 

At this point an important simplification occurs. Originally the path C„ is given 
by a complicated function z^(t) in four dimensions. However the integral over 
in (|65p is always from a point Z/j(t) on C„ to infinity and it does not depend on 
the particular form of C„ — it is the same as for the straight-line Polyakov loop. 
This is true, however, only for the Duki and not for Dikim- 

Another important observation is that the correlator D^^xp{u — v) is not a 
periodic function of (u4 — W4), in contrast to the fields F^,j{u) and F\a{v). This will 
be true also for the path integrals from any given point x to an arbitrary point 
y, and is a consequence of the vacuum average of the parallel transporter $(u,?;) 
present in D^yXp{u - v). 

Notice that the mixed contributions in (|64p . 



Jn" = \ j DiAki{u - v)da^i{u)d(TM{v), (67) 

is T- independent, so we shall neglect it in what follows. 

We finally turn to the spatial term, which will play a special role in what follows. 
Here the main contribution comes from the colour-magnetic correlator {u — v), 
which provides the area law for the (closed) spatial projection ^3 of the surface S'„. 
Correspondingly we shall write this term as (for ^3 3> , where Tg is the gluon 
correlation length, (x) ~ exp(— x/Tg)): 

(M/3(C„)) =exp(-(7,A3), (68) 

where the spatial string tension is defined in |T]) and A3 is the minimal area of the 
spacial projection of the surface Sn- Similarly to the case of the , one should 
start with white states of qq, gg or iq. As will be seen, the colour-magnetic vacuum 
{D^) acts in the centre-of-mass frame only in the states of the system with L 7^ 
and the resulting contribution of the does not separate into single-line terms, 
but it rather acts pairwise (triple- wise for iq) . Therefore one can account for colour- 
magnetic interaction in the higher (2-line or 3-line) terms. We shall neglect colour- 
magnetic interaction in the first approximation, as it is weak as compared to the 
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colour-electric one. As a result, the 4D dynamics in (|59p separates into 3D and ID 
ones, and it is possible to write 

G(")(s) = G^r\s)G3is) = J {Dzire-''--'"G3is) (69) 
and, similarly, 5*^"^(s) = 6*4" (3)5*3(3), with the free quark (gluon) factors 

One should notice at this point that G'"-' (s) and S'*^"-' (s) differ in the spin factors 
and in the colour group representation. Neglecting, in the lowest approximation, the 
spin factors, one obtains S'^"' from G^"-* by simply replac ing ad joint quantities by 
fundamental ones. Note that one can use Casimir scaling ^^1*^^ ! to write (Ts(adj) = 
|o-s(fund). 

To compute G^^"^ (s) one can take into account that does not depend on Z4 
and can be pulled out of the integral, while the remaining integral over Z4 can be 
taken exactly. In particular, splitting the proper-time interval Ne = s, ^ 00 one 
can write: 



rfP4^_,,,„^_ri.^ 1 (71) 
2tt 2V7rs 



As a result, one arrives at 



G^-\s) = -^e-^-^^Gr\s) (72) 



Substituting (|69p - ([7^ into (pT|) and (p^ . one arrives at the expressions for 
the pressure with the only unknown quantity being the fundamental Polyakov line 

. In the next chapter we discuss in detail the behaviour of the Polyakov lines at 
finite temperatures. 



3.3. Polyakov lines at finite temperatures 

An important result of the previous chapter is that colour-electric correlator Df 
yields factorised contribution, that is -^(12) = Z1Z2, where each individual factor 
Zi contains only the part of the common loop in the form of the Polyakov loop 
factor (with the singlet free energy in the exponent). The contribution of colour- 
magnetic fields does not factorize, however, this contribution can be considered as 
a correction, so it will be neglected. In addition, the average of the Polyakov line 
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operator was related to the colour-electric correlator jf' and to the static potential 

Vi(oo). In this chapter we discuss the behaviour of the Polyakov lines and static 

potentials in the Gaussian approximation to the QCD vacuum. 

At Tc > r > one has four Gaussian correlators D^{x), Df{x), D^{x), and 

D^{x), with the string tensions a^'^ given by H]). At T > the correlator 

(and, consequently, a^) vanishes, as was suggested in tSEI ^nd then was proved 

on the lattice El. The three remaining correlators are nonzero, moreover the spatial 

strin g te nsion as = cr^ grows with the temperature in the dimensionally reduced 
|oqi 

limit This fact e xplains also the growth with temperature of the Debye mass, 
mn — ^\f5's which is known from the lattice data Apart from this 

quantity, we shall not use below, in this section, the colour-magnetic correlators, 
since they do not produce static potentials between particles in the relative S'-wave. 
Therefore we are interested only in the colour-electric correlators D^{x) (inside the 
confined phase bounded by the curve Tc{p)) and -Of (x) (in the whole /x,T plane). 

Let us stress once more that, in our approach, only gauge- invariant states \n) 
are to be considered in the partition function ([63|) at T > 0. This is evident in 
the confined phase, since any coloured part of the given gauge-invariant system 
is connected to other parts by strings. In absence of colour-electric strings in the 
deconfined phase the necessity of using gauge-invariant amplitudes is less evident, 
except for world-lines in the spatial directions, where colour-magnetic confinement 
with a nonzero as is operating. Our use of gauge-invariant amplitudes, which fac- 
torise at large interparticle distances in the deconfined phase, leads to an explicit 
prediction for the EoS with modulus of phase factors, which approximately equal 
to modulus of Polyakov lines. 

Below we use, as in I 35 | 36 |^ gauge-invariant states, |n) at all temperatures and 
chemical potentials and describe the interparticle dynamics in terms of gauge- 
invariant quantities, like pair-wise or triple static potentials. The large-distance 
limit of these potentials yields one-particle characteristics — the self-energy parts 
of quarks, antiquarks, gluons, and so on. O ne can use those to study thermody- 
namics of QGP in the one-particle, or SLA I35 | 36| j|. jg rewarding, that the field 
correlator method is a natural instrument in describing this deconfined dynamics 
since, in absence of the , the correlator Df has the form of the full derivative and 
produces gauge-invariant one-particle pieces — self-energy parts — automatically 
(in addition to the interparticle interaction decreasing at large distances). 

Gauge-invariant states formed with the help of parallel transporters 

(Schwinger lines) create, as shown in Wilson loops W{C) for qq, qqq, or 
other quark systems. Latti ce data allow one to extract static interquark potentials 
in such multiquark systems I^^ESl. "When treating coloured systems, like (qq), taken 
as a part of a gauge-invariant system {qqqq in this case) , the pairs (qq) and {qq) are 
separated at a large distance and the potential V{qq, qq) is neglected. 

We start from the colour-singlet qq system and write contributions of the , 
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and Df to the static potentials at a nonzero T = 1/ (i^ 



Voir,T) = 2 di^il-iyT) {r-Od^D^{Ve+^), (73) 
Jq Jo 

r/3 



Jq Jq 
which contribute to the modulus of the Polyakov loops as 



+ ^^2), (74) 



Lfu.d = exp 



where Vi(T) = Vi{oo,T), Vq = Vo{r*,T), and r* is an average distance between 
the heavy-quark line and light antiquark (for > 0) or between the "hea yy-g luon 
line" and a gluon for Ladj ■ The Casimir scaling relation ([75]) predicted in is in 
good agreement with lattice data^^, as well as vanishing of the Lfund for T <Tc, 
Uf = and the strong decrease of Ladj for T < Tc- Indeed, for T < Tc and n/ = 
one has r* ^ oo and Vq ^ oo, explaining vanishing of the L^^j- For the l'"^} in 
this region one can take into account the kinetic energy of the gluon in the system 
consisting of an adjoint source plus a gluon (gluclump). This yiel ds an estimate 
LadjiT < Tc) — exp(— 77ig;p/r), where rrigip was computed in 1^ to be of order 
1 GeV. 

Notice that the Polyakov lines measured on the lattice are expressed through 
the (singlet) free energy of a QQ system at large distances i^i,^(oo, T), in the same 
way as in ((75|) . 

^fund = exp — . (76) 



In order to relate L^^'' and L^^^ one can use the standard representation for the 
free energy F^Q(r, T), 



exp I I = 2^ c„ exp I I , (77) 



where the summation covers all excited and bound states involving QQ, and 
(r, T) is the energy term of such a state n with the distance between the static 
charges Q and Q being equal to r. It is clear that L^J^ coincides with L^J^ when all 

states n, except for the ground state n = 0, are neglected. In this case V^'^ir^T) 
coincides with Vi{r,T) and, hence, with the F3,^(r,T). Note that Vi(r,T) in 
does not depend on T in the string limit of QCD, when the vacuum correlation 
length Tg tends to zero. 

In the general case all states n{QQ) contribute and therefore (c„ > 0) one has 
the following inequality: 

V^{r,T) > F^^ir^T) (78) 
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In order to define Vi and £fund properly one needs to separate their perturbative 
and nonperturbative parts and to renormalise Vi to g et rid of perimeter divergences. 
The separation in Df{x) can be seen at small 



Df{x) 



1 + 0{as In'^ x) TT^Gs 



247V, 



+ 



(79) 



and at large x, where (Z?f )"'' (x) behaves as 



{DfY'{x)^A,- 



-Ma\x\ 



where Mq is the lowest gluelump mass's], Mq w 1 GeV, Cp = (7V^ - l)/(27Vc) 
4/3 is the fundamental Casimir operator, and G2 is the gluon condensate, G2 
{as/n){F-,Fp. 

The corresponding separation of the Vi(r, T) is done in 1^ and reads: 

Vi{r, T) = Vl"\r, T) + Vl'^{r, T) + Vf^{a), 

where 

CpCts , 



Ai = 2CFasMoa: 



adj, 



(80) 



FP°"(r,r) = 



^(l + 0(rT)), 



(81) 
(82) 



is given by dH]), with the substitution Df (Df )"p, and 



div / 



O(rina) 



(83) 

and a is the 



Here mu = mjjiT) w 2^fal is the nonperturbative Debye mass 
lattice spacing (cut-off). 

The renormalisation procedure suggested in S2l amounts to discarding V^]^'^(a), 
and this is in agreement with the lattice renormalisation used in^Sl, where Fqq (f, T) 
was adjusted to the form V^°'^{r, T) at small r and T. Note, that V^^{r, T) - 0{r'^) 
in this region and the procedure indeed allows one to eliminate the constant term 
Vf^{a). 

Let us discuss now the contribution of the Vi to the interaction in QQ, QQ, 
and 3(5 systems. We start with the one-particle limit of the Vi{r, T) and the corre- 
sponding contribution to the if^j- According to the discussion above, one defines 
the renormalised Polyakov loop as in dill), ([81]), with Vi(T) = V^^{(X),T). We shall 
neglect the difference between l''^^ and ij^' (this difference becomes important at 
large T's, when > 1, while L^^^^ < 1). 



From (go]) one has (at T < 



so that, for A/q 



V^r(oo,T) = 
1 GeV, 

y7(oo,r,) 



Mi 



1 - 



T 



6asiMo)af 



l-e 



Mo 



0.5 GeV. 



(84) 



(85) 
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Fig. 1. Shown on the figure are curves of Ladj (blue dashed) and Lfund (red dashed) compared to 
the ones taken from In the T K Tc region the M{a.s = 0.195) = 0.982 GeV gluelump mass 
was used. In the deconfinement region the fit (1861 was used with Tc = 270 MeV for Lfund E'nd the 
Casimir scaled value for L^dj ■ 



A similar estimate one obtains from the lattice data which, at T ^ Tc, can be 
parameterize as 

F^Q(oo,r)« , T^Q((X3,T,)«0.5GeV. (86) 

1.35 (^) - 1 



Thus quarks (and antiquaries) have self-energy parts which, at T « Tc, are given 
by Kq{T) = Kq{T) = Vi{T)/2 « F^q(oo,T)/2 « 0.25 GeV. 

To illustrate the discussion of the Vb, Vi, and Tfund, ^adj we compare in Fig. [T] 
our theoret ical curves for the Tfund and Tadj (dashed curves) with the lattice data 
taken from 1^ (dots). The theoretical curves follow from (l75l) . with 14(00, T) = 
F^^(oo,T) taken from ([861) for T > Tc and Tadj = exp (-^) for T < Tc. From 
Fig. [H one can see a good agreement between our theoretical predictions and the 
lattice data. 

Similarly, for gluons, one has at T ~ Tc that Kg{T) = (9/4)Kq(T) « 0.56 GeV. 

Let us turn now to the r-dependence of the interaction. The perturbative part 
has a standard screened Coulomb behaviour (|82p . while the nonperturbative part 
vanishes at small r's: 

V"P(r,T) - const T^. (87) 

r— >0 



From ((111) and (HO]) one hasl- 

V^^ir^T) = V^'^i^^T) - ^Xi(Mor)Mor + O ) ^ n"'^(oo,T) + v{r,T). 



(88) 
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Fig. 2. The Casimir-scaled bare Polyakov loops for different representations D measured on 32'^ X 4 
lattices ESI. 



Therefore the nonperturbative interaction in the white system QQ changes from 
Vi"^(oo, T) Ri 0.5 GeV at large r's to zero at small r's. The same (multiplied by the 
factor 9/4) is true for the white gg system. 

We conclude this chapter by a discussion of the role of excited states in the 
definition of the F^q and of a possible violation of the Casimir scaling for the 
ifund and iadj- It is clear that in the Fqq for Uf ~ the only possible excited 
states consist of gluons: {Qg){Qg), {Qgg){Qgg), and so on. As it was shown inl3Sl^ 
the weakly bound states (Qg) are indeed supported by Vi (r, T) and, in neglect of 
the small binding energy, the total energy of these states is roughly given by the 
sum of the self-energy parts kq and Kg, 



Egg ~ ^^i(oo,T) + ^V^i(oo,r) « 0.8 GeV(r « T,). 



(89) 



This should be compared with the possible bound state of an adjoint static source 
G plus gluon which, in the weakly binding limit, can be written as 



EGg^2x -Fi(oo,T) « 1.1 GeV(T « TJ. 



(90) 



Notice that multiplicities of the states ([55]) and ([50]) are different, which leads to 
different predictions for corrections to the Fqq and Fqq. These corrections are 
not connected by the Casimir scaling, in contrast to the main (ground-state) term, 
for which V^q = Vi{oo,T) and Vqq = |Vi(oo,r). Therefore one expects some 
violation of the Casimir scaling by gluon-induced bound states in the Lfund and 

lool 

Ladj, though high precision lattice data ■ indicate a small role of such bound 
states. 
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Fig. 3. Pressure as function of temperature T. Shown on the left figure is a comparison 

of the analytical calculus l|95| lf dashed lines) with the lattice results (bold lines) HZl for the case 
Uj = 0,2,3. Shown on the right figure is the case ofrif = 2 + 1 . Green dashed line is the analytical 
calculation 1194 l l compared to the lattice one froml^SI 



3.4. Equation of state 

Here we follow the approach to study the QGP dynamics ^^^1^, where the main em- 
phasis was done on the vacuum fields, and the resulting modification of quark and 
gluon propagators was considered as the first and the basic step in the nonpertur- 
bative treatment of QGP — the SLA discussed before. As a result one obtains the 
nonperturbative EoS of QGP in the form of free quark and gluon terms multiplied 
by the vacuum induced factors. The latter are expressed via the only (nonconfin- 
ing) colour-electric correlator Df{x) and happened to be approximately equal to 
the absolute values of Polyakov loops Lf^nd, .^adj for quarks and gluons, respectively. 
Furthermore, due to the Casimir scaling property 3SJ£i^ these Polyakov loops are 
related to one another, as iadj — (-^fund)^^^- 

Moreover, the phase diagram was calculated in the SLA assuming that the phase 
transition is vacuum dominated, that is, a transition from the confining vacuum 
with vacuum energy density Econf — — ■^G2(conf) to the nonconfining vacuum 
with Edcc = —^02 (dec). The resulting phase curve Tc(/i) in depends on AG2 = 
G2 (conf ) — G2 (dec) and it was found to be in good agreement with lattice data 
for standard values of the G2(conf) and AG2 ~ O.5G2 — this will be discussed in 
detail in the next chapter. 

Thus the SLA is a reasonable starting point with no fittin g or model parameters, 
since Lfund can be computed analytically'^ or on the lattice^^, and AG2 « O.5G2 
is a fund ame ntal parameter of QCD. This picture of the QCD phase transit ion w as 
called in 1321 ^j^g Vacuum Dominance Model (VDM) originally proposed in™ in 
a simplified form (sometimes called the Evaporation Model). 

In this chapter we exploit the reduced pressure p = a-nd combine (|6H) , (|62p 
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Fig. 4. Pressure as function of temperature T. The case ol Uf = 2 + 1 (left figure) and n f 
(right figure) l|95p . Lattice results were taken from ESI 



and (|7ip .(f72 |) to obtain a simple form, which can be written as: 

pSLA AAT „ °° ( i\n+l 



n=l 

and 



^ _ P,f^ 2{Nl - 1) - 1 

^'Si - - —2 2^ ^^adj , (92) 



=1 

for quarks and gluons, respectively. Here (^^"^ is 

where an integral representation of the Bessel function K2 was used in the second 
equality. Then both sums (|9ip and ((92)) can be evaluated explicitly to yield: 

with a, = Vi(T)/2T, Ugi = fa,. 

Here we consider the case of /i = (generalisation to the case of nonzero chemi- 
cal potentials will be discussed in chapter [5]) and characteristic temperature region 
oiT K, Tc {Tc = 170 270 MeV) where quark masses do not affect the thermody- 
namical functions appreciably. This is due to the fast convergence of the sum over 
n at large n's ensured by the factors 1/n^ and L" (L < 1), while ipq^^ « 1 for rt ~ 1. 
Indeed, for — (^g"'' — 1, whereas for ruq — OAT one has (^q = 0.96, while 
fq'^^ = 0.03. Therefore one can neglect, with a good accuracy, the masses in (IMl) . 
arriving at: 

2nf r z^dz 8 r°° z^dz 



Zllf / z uz "/ z uz 

Pi^^ e-+a, + 1' Pa' Jo e^+-^'-l ^^^^ 



The results presented in (|95p are compared with the lattice pressure data, in 
Figs. [HI for Uf — 2 + 1 (left) and Uf — 3 (right figure). In Fig. [4] we show our 
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Fig. 5. Energy density as function of temperature T. The case nf =2 + 1 with rrt^ = 0.1ms 
and m„ ii = 0.2ms (green dashed curve) II97I I is compared to lattice data fromHSl(left fig.). The 
case rif = 2 + 1 with m^ = OAT, iris = T (red dashed curve) and rif = 3 with mq = OAT (blue 
dashed curve) II97I I are compared to lattice data from (right fig.). 



calculated curves for the cases n/ = 2 + 1 (left part) and Uf = 3 (right part), which 
are compared with lattice data from^^. 

As a further simplification, one can use, instead of the first terms of the 
expansion in (pi]) . (15^ . namely: 

127if 16 , , 

Pq = ^-^fund, Pg = ^^adj- (96) 

TT TT 

Other useful quantities to compare with the lattice data are the internal energy 
density and the "nonideality" of the QGP: 

<^ = T'4^(^) =e,+e,i (97) 



dT \ T 

and 



/(T)^^^t|, (98) 



respectively. 

Using ([94]) and ([96l) one has 

_ ^ 2 d dz 



(0) ^y-^y2j^ / ^ .gg^ 



and 



TT^ dJ TT^ dl 

The results of the calculations for are compared, in Fig. El El with three 
different sets of the lattice data:S3llHl^ andl^Il. In Fig.[7lwe demonstrate the /(T), 
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Fig. 6. Energy density as function of temperature T. The curve for = 3 with rrtu = 2 MeV, 
rrijj = 6 MeV, rus = 100 MeV (green dashed) I I97I I is compared to lattice data from EH 
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Fig. 7. "Nonideality" of QGP (e-3p)/T'*. Shown are the curves for (left fig.) = 3 with rriu = 2 
MeV, rrijj = 6 MeV, rUg = 100 MeV (green dashed line) compared to EH and (right fig.) for 
n/ =2 + 1 with rriu^d = 0.1ms and m„ = 0.2ms compared toHSI. Analytical calculations are 
done using II94I I and 1971 1. 



computed with the help of (jlOip and ([95)l . with lattice data for 2+1 flavours from 
1^ (left curve) and from SSI (right curve). 

It is instructive to estimate the contribution of the qq and gg interactions to the 
pressure. Writing the virial coefficient in the form = pf^(l + {pf^/T)B^{T) + 
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. . .), where Pj — Pq or Pgi in the SLA, with 

^'^^^ = ^ / " e^^^---^)/^) dV, J = fund, adj, (102) 

and taking the interaction terms C/fund and C/adj, at large T, as Uj{r^ T) = Tuj{rT), 
one obtains a corrected pressure in the form: 

P = Pf (l-c,) + p(°'(l-c,z) (103) 

where 

= p2dp(el"-^('')l - 1), cq^^ r p2dp(el«.w(p)l _ 1). (104) 

""Jo ""Jo 

Note that the and gg interactions in the singlet colour state is attractive, 
so that \Uj\ — —Uj. The dependence of the Uj on rT occurs at large T (in the 
dimensionally reduced regime) , when the dynamical dimensional quantities are the 
spatial string tension an = const • and the Debye mass mu{T) = 2^ an — 
const • T. 

Thus one expects that: i) the corrected pressure (|103p is smaller than the SLA 
prediction and ii) the large-T behaviour of the P{T) is below the Stefan-Boltzmann 
values (modulo logarithmic factors). Both features are clearly seen in Figs. [H [5] 



4. Beyond Single Line Approximation: Interacting quarks and 
gluons 

In the previous Sections one-particle contributions to the partition function were 
considered in the so-called SLA. It was stressed however that, at T > Tc, some 
interactions, like the residual colour-electric and colour- magnetic one, acts in white 
ensembles of quarks and gluons. Such interactions were neglected before, except 
for the contributions of the correlator Df which produced effectively single-line 
term Vi(oo). In ad dition, as well as Z)^, and produce bound systems of 
quarks and gluons | 13 | 42 | 46 | j^^^y affect strongly the dynamics of the QGP 

at T > Tc- Such bound states will be discussed below. To this end we start from 
a Green's function of the given quark(gluon) system and extract, in the same way 
as it was done before for the quark-antiquark meson at T = 0, the corresponding 
Hamiltonian i7, 

G(out|in) = (out|e-'f^/'^|in). (105) 

Our pm'pose in this Section is to discuss the properties of this resulting Hamilto- 
nian and possible consequences for the quark-gluon thermodynamics. Bound states 
in the qq, gg, and gq systems were studied in the framework of the FCM ^^-1 
and on the lattice — see and references in I 42 | 46 | Qj^g ^.^^^ consider two distinct 
dynamics: the colour-electric one, due to Df, and the colour-magnetic one, due to 
D" and Df . 
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iQQ)i 


iQQh 


iQQ)3 


{QQ)6 


iQg)3 


[Qgh 


(55)1 


i99)s 


a 





9/8 


1/2 


5/4 


1/2 


5/4 





9/8 


b 


1 


-1/8 


1/2 


-1/4 


9/8 


3/8 


9/4 


1 



Table 1. Parameters of the static potential of binary systems {QQ)d, {QQ)d, (Qg)D, and (gg)D 
in different colour representations D. 



4.1. Bound states at T > Tc due to colour- electric forces 

In the co' 
cussion) : 



In the colour-electric case the interaction can be written as (see | 4^ | 46 | Jqj. ^^iq dis- 



(r, T) = V, (r, T) = Vi (r, T) - (00, T) (106) 



where 



/■l/T rr 

Viir,T)= {l-vT)dv id^D^i^^y^Te)- (107) 

Jo Jo 

Notice that, in (|106p . we subtracted Vi(oo,T) which was already accounted in the 
SLA, in the form of the Polyakov lines — see ([55]) . 

Characteristic feature of ui(r, T) is that it is short-ranged, Vf^ff ~ 0.3 fm, and 
it can support bound 5-wave states in (cc)i, [gg)i, {cg)^, and {gg)s (the subscript 
stands for the representation of the colour group). The binding is weak, \e\ < 0.14 
GeV for T < 1.5rc. On the lattice, in addition to S'-wave states (cc)i, (66)1, and 
light (gg), the lowest P-wave state is claimed to exist (see for a review). 

It is clear that weakly bound states should dissociate fast in the dense QGP and 
therefore they hardly produce any significant effect on the production rate. They 
are important, however, for the kinetic coefficients. 

First, we generalise the static interaction (|107p to the case, when two colour 
objects, A and B, combine into a common colour state D 

yiAB,D)^^ r^^ = ^{C^^//W)(oo,r) + {Ca + Cb- CD)vl'^'^\r,T)} 



= aVl^^\oo,T) + bvl^^\r,T). (108) 

Here Cp = | is the Casimir operator for the fundamental charges, while Ca, Cb, 
and Cd are that for the representations A, B, and D, respectively. In Table [T] we 
illustrate ((108]). 



''The construction is similar to that found in 153] for the coefficient b of the vj^^^ (r, T) but differs 
in the presence of the first term proportional to a, since in EH the constant term was not taken 
into account. 
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For multicomponent systems one can similarly find, 



(109) 



^{Ca + Cb + Cc- Cd) E vI'^'^\u ~ r„T)}. 



In particular, 



^WQ)^ 1^^W0)(^^_^^. (110) 

i>j 

yiQQQ, ^ ^y('^'5)(oo,T)-i^y/^^)(r.-r„T) (111) 

yiQQQ) ^ ^ym}^^^T) + ly2vl''''\r.-r„T). (112) 

l>j 

It is easy to check that, at large distances, all systems consisting of ng quarks 
or antiquarks and Ug gluons tend to the constant limit, independent of _D, 

^^("qQ,",9)(|^^ - r,| ^ «3) = Equq + EgUg, (113) 



where 



EQ^\vl''^\^,n i?, = ^F/'3'5)(oo,r). (114) 



Since the nonperturbative part of the V^'^'^\r,T) behaves as O(r^) at r ^ 0, 
one obtains the lower bound on the nonperturbative part of l/("«'3^"99)(ry , T): 

a consequence, one can pre- 
dict the absence of bound states in some particular channels, for example, in {QQ)b,, 
{QQ)(ii {QQQ)w, and so on. 

As an application of the general relation (|87|) we show, in Fig. [H the static 
potentials y'-'^'^'^-'(ri,r2,r3,T) — V'^^\r^T) of three static fundamental quarks in 
three different representations D: singlet {Cd = 0), octet (Cd = 3) and decuplet, 
[Cd = 6). The quarks form a symmetric configuration of a equilateral triangle with 
the sides r. From Fig. [8] one can see that, indeed, all three potentials tend to the 
same limit (3/2)1//'^'^^ (r, T) at large r's, in agreement with (|113p . while deviations 
from this asymptotic at all distances are proportional to (1/2.-1/4,1/8) for the 
singlet, decuplet, and octet, respectively, as prescribed by (|110p - (|112p . This is in 
agreement with lattice calculations of free energies Fqqq (r, T) presented in 'SSI. 

Having constructed static potentials, we can now exploit the relativistic Hamil- 
tonian technique, developed in l ^3 | 21 | successful ly u sed for mesons, baryons, 
glueballs, and hybrids in the confinement phase (see ^^^^ for a review). This tech- 
nique does not take into account chiral degrees of freedom. Therefore, below we stick 
to the Hamiltonian technique of and consider heavy quarkonia and baryons. 
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Fig. 8. The three-quark potential (r, T) (in GeV) in different colour states for quarks 

placed at the vertices of an equilateral triangle with sides r versus r (in GeV~"^). 



Similarly to the case of the quark-antiquark meson discussed in detail before 
the bound-state problem for a generic multiquark(giuon) system can be formulated 
in the form (we introduce the einbein variables /i^ for quarks and gluons): 

H^^n = entlM, H = Ho + Hs + Hse, (115) 

where 

Ho= f^ + F/"«^'"«^)(r.,) (116) 

i=i "'^^ 

and Hs, Hse are the spin-dependent and self-energy parts of Hamiltonian defined 
in terms of field correlators. Then the mass of the system is given by 

E (^^ + yj+£n(m,-Mn«+nJ + a^/'^'^^(^,r)|, (117) 

after the minimisation performed with respect to all einbeins {/ii}. The stationary 
values {/i^^''} play the role of the constituent masses. In a search for bound 
quark-antiquark states was performed in the potential given by the sum of the 
nonperturbative colour-electric interaction and the screened Coulomb interaction 
(see dZll), dHOl), dMl) above and the discussion in 43), 

V^^{r,T) = F"P(r,T) + Vf{r,T). (118) 
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Fig. 9. Masses (first plot) and binding energies (second plot) of the cc colour-singlet bound states 
(in GeV) as functions of T/Tc- The curves are numbered in accordance with the data sets. 
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Fig. 10. Masses of bound states (in GeV) for the systems (cc)3, (cg)3, (gg)i, and (gg)s as functions 
of T/Tc 



For the gluelump mass Mq and the Debye mass — see ([M]) and — several 
data sets were used: in set I Mq = 1.0 GeV, to£i=0.69 GeV, in set II Mq = 1.0 
GeV, mD=0.2 GeV, in set III Mq 0.69 GeV, mD=0.69 GeV, in set IV Mq = 0.69 
GeV, mD=0.2 GeV, in set V Mq = 1.044 GeV, mD=0.2 GeV, and in set VI Mq = 
IMA GeV, m£, = 1.044 GeV. The results are presented at Figs. H [lOl From these 
figures one can see that, indeed, colour-electric interactions above the deconfinement 
temperature are strong enough to maintain bou nd s tates of quarks and gluons. 
Further details of the calculations can be found in SSl^ 
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4.2. Bound states at T > Tc due to colour-magnetic forces 

In this chapter we investigate bound states of quarks (generahsation to giuons is 
trivial) above the which app ear due to the residual nonperturbative colour- 
magnetic interactions. We follow ^12! 

We start from the Hamiltonian of a quark-antiquark meson (j40[) . Two particular 
cases are of most interest. The first such case corresponds to equal masses (mi = 
1712 — fn and therefore /ii = /i2 = A*), whereas in the other case one mass is assumed 
infinitely large (mi —^ oo, m2 = m and /ii — > oo, = m)- Then the Hamiltonian 
reads: 

H = - { — + M + / 4— +n+ '^2'' ^ ' 



(119) 

where ^ = 1 for the case of equal masses and ^ = 2 for the case of the heavy-light 
system. 

We take the extrema in and ry(/3) now, approximating ri{(3) by a uniform 
in P distribution. Then the Hamiltonian (|119p takes the form 

9^ + m^ 

and the spin-independent potential reads: 

1 „ \ _ „ , i..,2 „ y 



i^=7f^^^+A^)+^si(r), (120) 



Vsiir) = T]o<7Er + [ -qo ] crnr + -fiy , 770 = : , (121) 

V770 / 4 arcsin y 

with y being the solution of the transcendental equation 

^^y 



(JHT^ 4y \ \ (TH J J \r]o J anr 

The interested reader can find the detail s of a s imilar evaluation performed for the 
Hamiltonian (fTTOI) with as = cth = cf inl 25|55 | 

The remaining einbein fi is to be considered as the variational parameter to 
minimise the spectrum of the Hamiltonian (|120p . Obviously, the extremal value of 
H depends on quantum numbers and acquires two contributions: one coming from 
the current quark mass m and the other, purely dynamical, contribution coming 
from the mean value of the radial component of the momentum pr. 

It is instructive to pinpoint the difference in the potential (I12ip below and above 
the Tc. At small r's, the potential (|121|) turns to the centrifugal barrier L{L + 
l)/(^/xr^), whereas its large-r behaviour differs dramatically for the temperatures 
below and above the Tc. Indeed, the leading large-r contribution to the inter-quark 
potential corresponds to y ^ 1 and, for T < Tc, reads: 

Vconiir) = <JEr. (123) 

This is the linear confinement which is of a purely colou r-elect ric nature and which 
admits angular-momentum-dependent corrections (see^^^^^^. 
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In the deconfinement phase, at T > T^, the colour-electric part of the potential 
P2ip vanishes, the leading long-range term coming from the angular-momentum- 
dependent part of the interaction: 

Vsi{r) = —^ ^ + ---- 124 

Interestingly, in the deconfinement phase in absence of the confining potential, 
the spin-independent interaction becomes short-ranged decreasing as at large 
inter-quark separations. This feature means the full compensation of the centrifugal 
barrier which would naively behave as instead. The reason is obvious: at large 
inter-quark separations, the effective quark mass ^ is to be compared to the "mass" 
of the string ar. The bound-state problem solved in the potential (|123p gives a large 
value {pr) oc aE{r)-, so that, even for light (massless) quarks, their effective mass 
H appears quite large (/^ ^ m). On the contrary, for light quarks and in absence 
of the strong confining interaction (|123p . the values of fi are small (/x w m) and 
can be neglected as compared to the string contribution (Ti/r. This makes the spin- 
dependent terms in the effective inter-quark interaction important in this regime, 
as opposed to the confinement phase, where they give only small corrections to the 
bound states formed in the confining potential (|123p . 

Now we turn to the derivation of spin-dependent contributions to the inter- 
quark potential. The full set of spin-dependent interactions is given in (|43p . The 
leading spin-dependent term is the spin-orbit interaction (we omit contributions 
of ' which bring about short-range terms 0{r^'^)), 

SiLi _ 52^2 \ fldyo_ ^ 2dVi 
2/i2 J \r dr r dr 



where Vi (r) can be expressed through the colour-electric and colour-magnetic field 
correlators, 

i'-^^'-Tdrr dXD- ir, A), l!^ ^ -i T dr f dX (l - D^ir, A). 
r dr r Jo Jq r dr r Jg Jq \ r J 

(126) 

Only the Vi potential survives above the Tc, so that the resulting spin-dependent 
potential reads: 

where 

(K/3 - m') ^ f - e/2)^ = ^{^ + V'o) (- - yr^) , (128) 

and S = S\-^ for the light -light system, and S is the light-quark spin, for the 
heavy-light quarkonium. 

It follows from (|127p that above the deconfinement temperature, for the states 
with the total momentum J = i + S*, Sh > and the potential Vsoif) becomes 
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Fig. 11. The profile of the effective potential II130I I for m = 1 GeV (solid line), m = 2 GeV (dashed 
line), and m = 3 GeV (dotted line). 



attractive, with a possibility to maintain bound states. Furthermore, its slow de- 
crease as r — > cxD suggests that an infinite number of bound states exists, with the 
binding energies asymptotically approaching zero. Let us study these bound states 
in more detail. Hereafter in this chapter ~ and we use the notation a for the 
magnetic tension an- 

In view of an obvious similarity of the light-light and heavy-light cases (the dif- 
ference manifesting itself only in numerical coefficients), we investigate numerically 
only the light-light system, as a paradigmatic example. Furthermore, for r ^ Tg, 
the potential (|127[) does not depend on the form of the correlator since 

Finally, we neglect the perturbative part of the inter-quark interaction for it is 
screened to a large extent contributing to short-ranged forces only whereas the 
effect discussed in this work is essentially a long-ranged effect. 

Therefore we study the spectrum of bound states in the potential 

arcsiny U \ 2 



V{r)=[ CTr + /V 7 TTTTT, TTT^T^, (130) 

^' \ y arcsiny/ /ir(/i + 2(iy(/3 - 1/2)2)) ' ^ 

which is the sum of the spin-independent term p2ip and the spin-orbital term 
(fT27l) : y is the solution of (fT22l) with as = 0. In Fig. [IT] we plot the effective 
potential (I130|) for three values of the quark mass: m — IGeV, 2GeV, and SGeV. 

The resulting eigenenergy Snr-iifJ-) is added then to the free part of the Hamil- 
tonian ()119p . 

9 

711 

MnAp)^ — + ^ + e„„;(^), (131) 
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Parameter 


mb, GeV 


mc, GeV 


ms, GeV 


(T, GeV2 


Tg, fm 


Value 


4.8 


1.44 


0.22 


0.2 


0.2 



Table 2. The set of parameters used for the numerical evaluation. 





bb 


cc 


ss 


rir = 


-0.007 


-0.19 


-45 


1 — ' 


-5x10-'* 


-0.015 


-2.7 



Table 3. The binding energy E„^l = M„^l — 2m (in MeV) for the ground state and for the first 
radial excitation in the potential (11301 1 with L = 1 for the bb, cc, and ss quarkonia. 
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Fig. 12. The binding energy of the quark-antiquark system versus the mass of the quark for L = I 
and Ur = (first plot) and = 1 (second plot). 



and this sum is minimised with respect to the einbein fi, 
5M„,/(/i) 



= 0, M„,.i = MnA^^o)■ (132) 

In Table [5] we present the set of parameters used in our numerical calculations, 
whereas in Table |3] we give the results for the binding energy for the bb, cc, and 
ss quarkonia above the Tc for I = 1 and n,. = 0, 1. We ensure therefore that for 
L ^ the potential (|130p does support bound states. The binding energy is small 
{\En^L\ T for the 6 and c quarks and I-Ekj-lI < T for s quarks) so these bound 
states can dissociate easily. 

Let us discuss the problem of binding of light quarks. The effective poten- 
tial (|130p admits different forms at different inter-quark separations, depending on 
which contribution, of the quark mass term ^ or of the "string mass" 2(i^(/3— 1/2)^)), 
gives the dominating contribution, that is for (x ^ ar and /i ^ or. If large 
distances contribute most to the bound state formation (the latter case), then 
V{r) — 0{L{L + l)/(oT^)) -I- 0{L/{fir'^)), where the first term comes from the 
spin-independent interaction (see (|126|) 1 and the other stems from the spin-orbit 
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potential. The dependence of the binding energy on /i is expected then to be rather 
moderate, approximately as 

On the contrary, in the former case with the string dynamics giving a correction 
to the quark mass term, the potential (|130[) can be approximated as 

L{L + 1) al 
V{r) w 133 

that is by the sum of the centrifugal barrier and the attractive Coulomb-like po- 
tential with the effective coupling 

cL 

ttoff = 134 
/•« 

The corresponding eigenenergy can be found in any textbook in Quantum Mechan- 
ics and gives a stronger dependence on /z, 

0-2^2 

e{fi) cx —^alff cx —. (135) 

M 

Let us consider the states with L = 1 and rir = 0- We follow now the procedure 
described in detail before, that is we solve the full problem numerically, and find 
the dependence of the eigenenergy e on the einbcin fi to be 

« (136) 

^^ 

Comparing this to (|135p we find a good agreement, with the small deviation in 
the power resulting from the proper string dynamics. We conclude therefore that 
the dynamics of the system develops at the inter-quark separations Tg ^ r < m/a 
(since, for the set of parameters given in Tabled m > uTg then there is room for 
such separations). 

Let us discuss now the procedure of minimisation of the spectrum (|13ip in /i. 
First of all, let us notice that, in the einbein field formalism, the calculation of the 
spectrum naively looks like a nonrelativistic calculation due to the "nonrelativistic" 
form of the kinetic energy in the Hamiltonian with the einbein field ^ introduced. In 
the meantime, the full relativistic form of the quark kinetic energy is readily restored 
as /X takes its extremal value and hence this is the procedure of taking extremum 
in fjL in the masses (|131|) to sum up an infinite series of relativistic corrections and 
thus to restore the relativistic spectrum. For example, the relativistic ground state 
eigenenergy Eq = rriy/l — (Za)'^ of the one-body Dirac equation with the Coulomb 
potential —Za/r can be reproduced exactly with the help of the einbein technique. 
Finally, one can visualise the form of fi considering the effective Dirac equation 
for the light quark in the field of the static antiquark source. When written in 
the form of a second-order differential equation, it contains the spin-orbit term 
of the same form as given in (|125p but with fi replaced by the combination e -I- 
m + U ~ V, where U and V are scalar and vector potentials, respectively. For 
light (massless) quarks this combination takes drastically different values below 
and above the deconfinement temperature. Indeed, in the confining phase of QCD, 
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when spontaneous chiral symmetry breaking leads to a strong effective, dynamically 
generated scalar potential U, this effective "/i" is large. On the contrary, above the 
Tc, when U is small "/z" is also small (it can be even negative since the eigenenergy 
e may have any sign). 

It has been found numerically that the extremum in /i for Af„^i(/i), (jl3ip . exists 
for m exceeding the value of approximately 0.22GeV (for the given a = 0.2GeV), 
and no extremum exists for smaller values of the quark mass (see Fig. [T^ for the 
dependence of the binding energy En^L on the quark mass). This property of the 
bound state spectrum can be easily understood using the analogy with the bound 
state problem for the Dirac equation with the potential in the form of a deep 
square well or Coulomb potential discussed above. For example, for the Coulomb 
potential, a problem appears as the coupling exceeds unity — the well-known 
problem of Z > 137. From (|134p we easily find this critical phenomenon to happen 
at m « /X < y/a « 0.4GeV. This estimate is in good agreement with the result of 
our direct numerical calculations quoted above. 

Physically this situation means that many quark-antiquark and/or gluon pairs 
are formed and finally stabilise the vacuum. Formally the problem is not anymore 
a two-body problem, but rather many-body, so that many-body techniques are 
to be applied. For example, in electrodynamics with Z>^137, one can derive the 
resulting self-consistent field of the Thomas-Fermi type A similar situation can 
be expected in the deconfinement phase of QCD. In absence of the linear potential, 
the einbein /x (playing the role of the effective quark mass) is not anymore bounded 
from below by the values of order ./oe — 0.4GeV coming from the binding energy in 
the linearly rising potential. To see the onset of this phenomenon in the framework 
of our two-body (one-body for the heavy-light case) Hamiltonian, one should take 
into account the negative-energy part of the spectrum, when the full matrix form of 
the Hamiltonian is considered Indeed, the matrix structure of the Hamiltonian 
occurs in the path-integral formalism from the two-fold time-forward/backward 
motion described by the positive/negative values of ii. Off-diagonal terms in the 
matrix Hamiltonian produce the turning points in the particle trajectory and result 
in Z-graphs. 

Notice that the same is true for the glueballs and gluelumps since in this case 
equations are the same as for light-light and heavy-light quarkonia, respectively, 
but with the quark spin replaced by the gluon spin and an by j(Th- 

It is important to notice that a separated quark-antiquark pair was considered 
in this chapter. In reality such quark-antiquark pairs are to be considered in the 
medium formed by other quarks and gluons, that is as a part of the SQGP. As a 
measure of the interaction in SQGP one can consider the ratio of the mean potential 
energy to the mean kinetic energy of the particles in the plasma, F = {V) / {K). It is 
easy to estimate that (A') ~ T and (V) ~ an/T. This gives F = an/T'^ and so this 
parameter is large for quarks and it is several times larger for gluons. Therefore, 
SQGP is a strongly interacting medium which looks like a liquid, rather then as a 
gas. With the growth of the temperature the medium becomes more dense, and the 
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mean distance between particles decreases. As this distance becomes comparable 
to the radius of the bound states discussed in this chapter, the latter will dissociate 
because of the screening effects. In other words, the hot medium plays the role of a 
natural cut-off for the effect of bound pair creation discussed above. Notice however 
that, for the quark masses around 0.2GeV, the radius of the bound state is of the 
order of one fm and it is expected to decrease further with the decrease of the quark 
mass, even if the pair creation process is properly taken into account. This means 
that indeed there is room for such bound states for the temperatures above the 
Tc- Breakup, with the growth of the temperature, of such high-Z states for quarks, 
and especially for gluons which possess more degrees of freedom than quarks, may 
affect such characteristics of the plasma as its free energy and it entropy (for a 
recent paper of explaining the near-Tc behaviour of these characteristics seei^^. 

Concluding this section one can say that colour-magnetic (spin-dependent) in- 
teraction acting on light quark or gluonic systems enforces nonperturbative creation 
of light qq and gg pairs. 



5. Generalisation to nonzero baryon densities 

In this section we will extend the method to the case of nonzero /i and nonzero 
interaction with vacuum {Vi ^ 0) and to find Tdp), for /i > 0. The simple pic- 
ture of VD dynamics with the only input Lf^j^n^lT), taken from lattice or analytic 
calculation suggested in^SSl^ which was discussed below as the Vacuum Dominance 
Model (VDM), is adopted below and is shown to produce surprisingly reasonable 
results, being in good agreement with available lattice data for nonzero /i, where 
these data are reliable. 



5.1. Nonperturbative EoS for fi > 

The main idea of the VDM is that the most important part of quark and gluon 
dynamics in the strongly interacting plasma is the interaction of each individual 
quark or gluon with vacuum fields. This interaction is derived from field correlators 
and is rigorously proved to be embodied in factors, which happen to coincide with 
the modulus of Polyakov loofjf], 

ifu„d = exp( — 1, Ladj=expl-- — I, (137) 

where Vi{T) = Vi{oo,T), Vd = Voir* ,T) and Viir,T),VD found in 13^ to be 
(/3 EE l/T) the same as for ^ = case (see Eqs. ([731) , (HH) ) but with , Df in 
principle depending on /i. 



'^We neglect in this approximation the difference between Lfund expressed via Vi(oo, T) and ifJJ'j^j 
where the role of Vi(T) is played ESEZIj^y singlet QQ free energy F?, = (oo, T). The latter quantity 
contains all excited states, so that Vi{T) > F3j.{oo,T). 



2, 2009 16:39 WSPC/INSTRUCTION FILE QGP 



40 A.V. Nefediev, Yu.A. Simonov, M.A. Trusov 

In Vd, Ea. p37p r* is the average size of the heavy- Hght Qq or its adjoint equiv- 
alent system; for T > Tc one has ^ Vd = 0. For T < Tc one has r^^j^^ = oo 
for Uf — 0, yielding £fund = 0, however r*^^ w 0.4 fm for any Uf and gives nonzero 
LadjiT < Tc). The form (|137p is in good agreement with lattice data'^ and also 
explains why Lfund is a good order parameter (however approximate for n/ ^ 0). 
Note, that only NP parts Df,D^ enter in ([72]) , dZ!]) , seeSlfor discussion of sepa- 
ration of these parts; renormalisation procedure is discussed also in 

In the lowest NP approximation one neglects pair, triple etc. interactions be- 
tween quarks and gluons (which are important for Tc < T <■ 1. 2Tc where density 
is low and screening by medium is not yet o peratin g, see l ^^ l ^^ ) and derives the 
following EoS (this approximation is called inl^^^l^the Single Line Approximation 
(SLA)) 

2{Nl - 1) - 1 

P9I = -J^- - ^2 2^ ;^^adj (139) 



=1 



with 



vS"'m-^''■=(^)-l-I(=?)^- ("») 

In inSD, it was assumed that T < i = 1 GeV, where A is the vacuum 

correlation length, e.g. d[^\x) ^ e"'^'/'*', hence powers of L", see^Sfor details. 

With few percent accuracy one can replace the sum in (jl39p by the first term, 
n—1, and this form will be used below for pgi, while for pq this replacement is not 
valid for large ^, and one can use instead the form equivalent to (|138p . 




(141) 



where v — mq/T and 



$,(a) = / -^^^—^ . (142) 

Eqs. (|14H) . (|139p define Pq,Pgi for all T, /i and mq, which is the current (pole) 
quark mass at the scale of the order of T. 

To draw pq , pgi and p = pq + pgi as functions of T, /i one needs explicit form of 
Vi (T) . This was obtained analytically and discussed inHH, another form was found 
from Df{x) measured on the lattice in|S2l and is given in SSI. 

Also from lattice correlator studies '^'^^J Vi{T = Tc) is (with ^ 10% accuracy) 
0.5 GeV and is decreasing with the growth of T (cf. Fig. 2 of 221 and Fig. 1 of SSI). 
This behaviour is similar to that found repeatedly on the lattice direct measurement 
of see e.g. Fig. 2 of where Fii^ = Vi{T) is given for n/ = 0,2,3. In what 
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follows we shall exploit the latter curves parameterizing them for T > Tc and all 
n/ as 

ViiT)^ "'^^^"^ ' Vi{Te)^0.5GeY. (143) 
1.35 f^j - 1 

For fi > one can expect a /i-dependence of Vi , however it should be weak for values 
of fi much smaller than the scale of change of vacuum fields. The latter scale can be 
identified with the dilaton mass Wd, which is of the order of the lowest glueball mass, 
i.e. ~ 1.5 GeV = m^. Hence one can expect, that Vi in the lowest app roxi mation 
does not depend on fi. This is supported by the lattice measurements in^^, where 
for Y' — 1-5 and ^ = 0.8 the values of F^'^ are almost indistinguishable from the 
case of ^ = 0. 

To give an illustration of the resulting EoS we draw in Fig. [T3] the pressure p 
for the cases n — 0, 0.2, 0.4 GeV and Uf — 2. One can see a reasonable behaviour 
similar to the lattice data, see for a review and discussion. 

5.2. Phase transition for nonzero /j, 

Here we extend the VD mechanism suggested in'^to the case of nonzero fi and Vi. 
We assume as was said in Introduction that the phase transition occurs from the full 
confining vacuum with all correlators , Df , , present to the deconfined 
vacuum where vanishes. The basics of our physical picture is that all fields (and 
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correlators) do not change for /x, T changing in a wide interval, unless u, T become 
comparable to the dilaton mass « M(glueball 0+) « 1.5 GeM^- Therefore 
correlators and , are almost constant till T = and at T > Tc a new 
vacuum phase with = = is realised, whi ch yi elds lower thermodynamic 
potential (higher pressure). Lattice measurements IBl^ support this picture. The 
crucial step is that one should take into account in the free energy F of the system 
also the free energy of the vacuum, i.e. vacuum energy density e^ac = j^i^^f^ = 
^e"'^ ii-^flv)'^) ~ ~'^~~l32~^^2: which can be estimated via the standard gluonic 
condensate IS21g2 = ^{{Fp^) w 0.012 GeV. 

Hence for the pressure P — —F one can write in the phase / (confined) 

Pi = \e.ac\+Xi{T) (144) 

where x(^) is the hadronic gas pressure, starting with pions, Xpion — 

Ir^T^.ln the 

deconfined phase one can write 



^dec I 



{Pgl+Pq)T^ (145) 



where le^^cl the vacuum energy dei isity in the deconfined phase, which is mostly 
(apart from -Df (0) ~ 0.2Z?^(0), see ^SEllj colourmagnetic energy density and by 
the same reasoning as before we take it as for T = 0, i.e. le^^cl — 0.5|ei,ac|- 
Equalizing Pj and P// at T = Tc(/i) one obtains the equation for 

T.(,)^r^'--'+^(^)V'' (146) 



Pgl + Pq 

where A|e„ac| = \£vac\ - \£tTc\ ~ AG2; AG2 ~ ^G2; Pgi and Pq are given 

in (|139p . (|141[) respectively and depend on both Tc and ^. 

In this letter we shall consider the simplest case when the contribution of 
hadronic gas Xi(r) can be neglected in the first approximation. Indeed, pionic 
gas yields only ~ 7% correction to the numerator of (|146[) at T « Tc, and from 

one concludes that x{Tc) ^ 0.5Tc , which yields a < 10% increase of Tc for 
G2 - 0.01 GeV. 

From the expression for Tc(//) (I146|) one can find limiting behaviour of T'c(a* ~* 0) 
and /Zc(T' — > 0). For the first one can use for pq and pg ()138|) and ()139p and expand 
r.h.s. of (jl46p in ratio Pg/pq with the result. 



Tc = t(") 1 + + ^li^ (147) 




"^Note that G2 does not depend on fi,nf in the leading order of the 1/Nc expansion and one 
expects a growth of the magnetic part of G2 at T > T^im.red ~ 2Tc, G™"^ Ri 0{T''), in the 
regime of dimensional reduction. 



2, 2009 16:39 WSPC/INSTRUCTION FILE QGP 



Deconfinement and quark- gluon plasma 43 



where the last term yields a 3% correction, and T^*^^ = ( ^^^ ■ Solving 



([TTf)) for Tc one has 

with K = iVi(Tc). From p46p . (|147p one can compute expansion rc(M) in powers 
of /i, 

Tc(mb) = r,(0) (^1 - , MB = 3m. 

C= =0.0110(3) for nf = 2,3,4. 

144VTT#f ^ ^ ^ ' ' 

One can see, that C practically does not depe nd on r if and is in the same ball- 
park as the values found by lattice calculations, see^^^^^for reviews and references. 

Another end point of the phase curve, ^c{T — > 0), is found from (|146p when one 
takes into account asymptotic <i>o(a ^ oo) = ^ + ^^a? + + which yields (for 
small ^) 

The resulting curve Tc(/x) according to Ea. p46p with xi = is given on Fig.fTH 
(left side) for AG2 = 0.00341GeV'*, n/ = 2,3 and rn„ = 0. To compare, we have 
shown on the same figure (right side) the lattice data^^ obtained in the reweighting 
technique. 




Fig. 14. The phase transition curve Tc(/x) from EQ. ljl46|l (in GeV) as function of quark chemical 
potential (in GeV) for rif = 2 (upper curve ) and rif = Z (lower curve) and AG2 = 0.0034 
GeV* (left side) in comparison with lattice results EZI (right side) 
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Table 4. The values of Tc{fJ. = 0) and fJ,c{T = 0) computed using 11481 1 and 11491 1 for several values 
of AG2 and = 0,2,3. 



AG2/(0.01 GeV"*) 


0.191 


0.341 


0.57 


1 


Tc(GeV) n/ = 


0.246 


0.273 


0.298 


0.328 


Tc(GeV) n/ = 2 


0.168 


0.19 


0.21 


0.236 


Tc(GeV) n/ = 3 


0.154 


0.172 


0.191 


0.214 


/ic(GeV) uf = 2 


0.576 


0.626 


0.68 


0.742 


/ic(GeV) = 3 


0.539 


0.581 


0.629 


0.686 



5.3. Discussion of results 

The prediction of Tc{^) depends only on two numbers: 1) the value of gluonic 
condensate AG2 and 2) the value of Vi{Tc) = 0.5 GeV taken from lattice data^Sl 
(and quantitatively close to the value from the analytic form' — J). 

We take G2 in the limits 0.004 GeV < G2 < 0.015 GeV*, the value G2 = 
0.008GeV4 (AG2 = 0.0034 GeV) being in agreement with lattice data of Tc{{)), 
for n/ = 0,2,3, see Table H 

Note that Tc at = in Table |4] is obtained not from (|148|) , but directly from 
([T461) with n/ = 0,x(T) = 0. 

The curve in Fig. [14] has the expected form, which agrees with the curve, ob- 
tained in^SZl \yy the reweighting technique and agrees for /i < 300 MeV with that, 
obtained by the density of state method and by the imaginary /j, method 

An analysis of the integral (|142p for v = Q reveals that it has a mild singular 
point at ^ging = ^±z7rT, which may show up in derivatives in ^. At T = 0, ^sing ~ 

= 0.25 GeV and is close to the point where one expects irregularities on the 
phase diagram ISSESI 

The limit of small ji is given in (|148p . Taking Vi{Tc) ~ 0.5 GeV as follows from 
lattice and analytic estimates, one obtains with ~ 3% accuracy the values of 
given on Fig. [Hfor n/ = 2, 3 and AG2 = 0.0034 GeV^. 
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The values of /ic, from (jl49p are given in Tableland are in agreement with the 
curves shown in Fig. [Ml Note that x(r = 0) = and (|149p holds also in the case, 
when hadron (and baryon) gas is taken into account. 

At this point one should stress that our calculation in VDM of p{T) does not 
contain model parameters and the only approximation is the neglect of interparticle 
interaction as compared to the interaction of each one with the vacuum (apart from 
neglect of x(T')). Fig. [T3l demonstrates that this approximation is reasonably good 
and one expects some 10—15% accuracy in prediction of Tc(/x). Note that in VDM 
the phase transition is of the first order, which is supported fo r n f = by lattice 
data, see e.g.l^however for n/ = 2, 3 lattice results disagree, see^^^lZllfoi- a possible 
preference of the first order transition. One however should have in mind that the 
final conclusion for lattice data depends on input quark masses and continuum 
limit. 

The "weakening" of the phase transition for n/ > and nonzero quark masses 
is explained in our approach by the flattening of the curve P{T) at T « Tc when 
hadronic gas x{T) is taken into account, since x(2^) for nj — Q la much smaller 
than for Uf — 2. 

Finally, as seen from our expressions for Pq,Pg (|138p . (|14ip . our EoS is indepen- 
dent of Z{Nc) factors and Z{Nc) symmetry is irrelevant for the NP dynamics in our 
approach. This result is a consequence of more general property - the gauge invari- 
ance of the partition function for all T, /i, which requires that only closed Wilson 
loops appear in the resulting expressions, yielding finally only absolute value of the 
Polyakov loop, or in other words, is expressed via only singlet free energy of quark 
and antiquark F^^J. « Vi. 

6. Concluding discussions 

6.1. Comparison to other approaches 

Our resulting EoS, in particular P{T, /i) in the leading (single line) approximation 
has the form of ideal gas factors multiplied by moduli of Polyakov Lines (PL) . In 
the next orders also interaction between quark and gluon lines is included, which 
leads to the decreasing of pressure, so the final results are in better agreement with 
lattice data, see Figs. 3-6. 

The interaction between quarks and gluons has both colourelectric and colour- 
magnetic components and was discussed in section 4. 

Thus it is clear that PL play a very important dynamical role in the FCM. This 
role is in creation of T-dependent self-energy term for each individual quark or 
gluon, which suppresses their contribution to the pressure. One can say therefore, 
that quarks and gluons acquire vacuum induced effective masses. 

From this point of view it is reasonable to discuss other approaches in the 
literature, where PL play important role. 

One should start with the formalism of effective action for Polyakov lines, ex- 
isting for almost 30 years see also for later development. 
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Here PL and Z{N) symmetry are mostly used conceptually as an order param- 
eter for the deconfinement transition, and the effective potential was calculated for 
purely gluon and quark-gluon system to study potential extrema in terms of PL 
moduli and phases. From the discussion in sections 1-4 it is clear however, that 
only moduli of PL enter EoS in our approach and Z(N) symmetry is irrelevant for 
dynamical properties of QGP. Moreover, quarks violate Z(N) symmetry and the 
latter has no connection with the deconfinement in the unquenched case. Instead, 
moduli of PL are crucially important as dynamical input in establishing EoS and 
the QCD phase diagram at nonzero fi. 

Indeed, it is the modulus of L{T), which defines in Eq. (|146p the transition 
temperature 71, (^) in the SLA. ^^^^ 

There is another approximation l ^6 | 77 | ^^le theory of QGP, which is closer to 
ours, since it also uses PL as a kind of dynamical input. It is called PNJL and 
is based on the Nambu-Jona-Lasinio Lagrangian with quark operators augmented 
by PL. Here one obtains EoS and effective quark and gluon masses in qualitative 
agreement with lattice data. However, since NJL Lagrangian is not confining, the 
detailed picture of deconfinement may be quite different. 

The action density of the PNJL model can be written as (see Rossner et al in 

for a review and additional references) 



where N is quartic in ijj as in NJL model, and ^4 is constant field entering PL, 
which therefore become complex in general. This fact and the associated Z{Nc) 
symmetry differs from our results and our EoS, where only modulus of PL enters. 
Nevertheless one can obtain in this appr oach a reasonable agreement for the first 
moments of pressure with respect to ^/T^^. 

At this point it is useful to compare our results with the phenomenological 
models of QGP, where qua rk and gluon masses have been introduced as input and 
fittin g param eters ,see e.g. I ^^ l ^^ l. 

In I '^^ 1 ^9 1 the quasiparticle propagators are introduced with real and imaginary 
parts of self-energy (effective mass) terms, which can be compared to our PL se\i- 

energies,, £1 = \ciVi{oo,T), Li = expi^-^^^'W^^i) ^ ^ = 1^9^ Ci = l,|. 

It is interesting that the resulting effective masses for quarks and gluons ap- 
proximately satisfy Casimir scaling as in our CiVi{oo, T), and have the same order 
of magnitude as in our calculations. A more detailed comparison will be subject of 
future publications. 

We now turn to the important question of the density induced phase transition, 
which may be important for several planned ion machines and possible existence 
of quark cores in neutron stars. As was discussed above in section 5, FCM predicts 
deconfinement phase transition at large fiq = ^Xcr ~ 0.6 GeV, provided vacuum 
properties and in particular gluon condensate and PL are not affected by jig. It is 
important to study what happens in the region /ijv l£ fJ-q l£ fJ'cr, where is the 
chemical potential at normal nuclear density, fiN ~ 0.3 GeV. 



S = ^|J+{^r - iaS/ + 747710 - ^4)7/- + N{^|J, 



(150) 
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It is clear that nonzero /i„ can influence baryon structure, and some change in 
the nucleon mass and radius was studied before . However much more drastic 
change due to fiq appears in the interquark forces, as was found iniSli. Here it was 
reaUzed that hght quark interaction and confinement is obtained self-consistently 
from nonlinear equations, yielding linear scalar confinement for vanishing but 
a deeper well for r < ^ appears for nonzero /i^, which makes baryon lighter and 
more compact. Moreover, 3nq systems may become more stable for n = 1, 2, ... due 
to larger number of pairing interquark forces. More detailed quantitative analysis is 
needed for the final positive answer, however it seems feasible , that at growing fiq 
the nuclear matter becomes a dense medium of smaller size 3q with an admixture 
of 6q,9q,... bags, the latter is growing in percentage with /z^, unless quarks find 
themselves finally in the deconfincd phase. In this scenario the role of multiquark 
bags and, hence, of cumulant mechanism is growing with density, which can be 
measured experimentally in future ion- ion experiments. 

Till now nothing was said about qq pairing and possible quark superconductivity, 
widely discussed in literature — see for reviews. The np vacuum changes this 
picture in several respects. The first one was discussed ini^, where it was stressed, 
that the minimum of thermodynamic potential Qq (with the gap in the scalar 
diquark sector Aq — (0.1 0.15) GeV is of the order of 10"'^ GeV* or less, and this 
value is of the order oinp vacuum density \s\, so the realization of the perturbative 
diquark condensation crucially depends on this np value and is questionable when 
l^ol < \e\. 

Another np feature of the QCD vacuum, which may prevent standard diquark 
scenario, is the qq interaction in the diquark sector, discussed in section 4. As can 
be seen in Eq. (jlOSp interaction in the qq sector can be written as 

Vqq^adVi{^,T) + bdVi{r,T), (151) 

where 03 = ^, ag = |, 63 = i, fog = —|; to be compared with qq system, where 
ai = 0, 61 = 1 and Vi{oo, T^) w 0.5 GeV. 

Here Vi{r,T) is the sum of perturbative and np interactions, Vi{r,T) = 
Vi(oo, T) + v{r, T), and Vnp{r, T) is negative everywhere. 

As a result one has additional (constant) repulsion in the qq system as compared 
to qq, of the order of iVi(oo,Tc) « 0.25 GeV, which is larger than the assumed 
gap Aq. Therefore one can expect that pairing is destroyed and one has only white 
correlations in the qq and qqq systems which can appear as bound states, e.g. in 
the S -wave cc system, which is supported by lattice data . 

Thus we are coming again to the " quark-bag scenario" discussed above, where 
nucleons in nuclear matter become tighter and lighter with increasing density and 
are accompanied with growing number of 6q, 9g, . . . quark bags finally occupying 
the whole volume. 

This scenario possibly implies a more hard EoS of nuclear matter, than usually 
exploited and can lead to higher masses of neutron stars with quark matter cores. 
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see in this respect the discussion in l ^5|86 | 

Finally, we discuss in interesting idea of Dreniin and collaborators, about the 
possibility of Cherenkov emission of gluons in QGP due to the effective creation 
of nonzero dielectric permittivity e in a dense medium .As was shown there, 
this Cherenkov emission can possibly explain the double-humped structure of the 
away-side jet at RHIC. We can add here from the point of view of FCM, that e is 
sensitive not only to the medium effects, but can also get contributions from the 
change of QCD vacuum properties, which happens at T > Tc- This topic is now 
under investigation. 

6.2. Mysteries of deconfinement 

The theory discussed in this review was able to address most important questions 
related to the QGP dynamics, namely: 

(1) Why does deconfinement take place at some temperature and density? The the- 
ory based on the FCM shows that thermodynamic potential (pressure) grows 
faster with the temperature T in the deconfined phase, as compared to the con- 
fined phase, so that the (almost) constant difference in vacuum energy density 
between the phases is overcome at some T ~ Tc. 

(2) What is the main dynamical difference between the phases? It is shown in the 
review that, in the confined phase, both colourelectric (CE) and colourmagnetic 
(CM) confinement coexist while, in the deconfined phase, only the latter is 
retained In the meantime, CE forces due to the nonconfining correlator Df 
also survive in the deconfined phase and they can support bound states at not 
too large T's, T < l.ST^. 

(3) What are manifestations of the CM confinement? As was discussed in this 
review, the CM confinement ensures strong interaction in the quark-gluon sys- 
tems and can support weakly bound states for nonzero angular momentum L. 
Moreover, most part of the spin-dependent forces is due to the CM confining 
correlator and thus it survives the deconfinement transition 1^3^ xhe CM 
confinement solves the old Linde paradox since the 3D perturbative diagrams 
acquire Wilson loops of the confining background field. The latter converge 
at large distances as {W) ^ exp(— cj^f x area) and make all integrals conver- 
gent However, this is probably only a small part of the story. Indeed, the 
very notion of confinement implies that, when moving with some velocity, any 
coloured subsystem is not free but is rather connected, by a CM string, to its 
colour partner, which makes the entire system white. As a consequence, a gluon 
or a quark acquires a nonperturbative CM mass, which is velocity-dependent. 
This effect was studied at the example of the Debye mass, which was calculated 
in and was shown to reproduce exactly the lattice data (in contrast to 
a pertur bati ve treatment, which is not gauge invariant and is about 50% off 
the data^^^. It is clear, that this CM dynamics is of a general character and 
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may be related to many features ascribed to QGP, for example to jet quench- 
ing, collective phenomena, and so on. This important topic requires a further 
careful investigation. 

(4) What is the connection between deconfinement and Chiral Symmetry (CS) 
restoration? On the theoretical side, it was found in the FCM that the same 
correlator , which ensures the CE confinement, enters the kernel of the 
nonlinear equation, th e latter demonstrating Spontaneous Breaking of Chiral 
Symmetry (SBCS) Thus confinement entails SBCS, and both should dis- 
appear at the same temperature. This behaviour is observed in most of lattice 
simulations (see, for example, ^^^1) However, in some papers (for example, in 

a difference AT ^ 20 MeV is found between the temperatures of the de- 
confinement and CS restoration, though it is not clear whether it is a physical 
effect or just a result of lattice approximations. One might inquire at this point 
of the role played by the CM confinement. This problem can be connected with 
the value of the chiral condensate {qq), as a function of temperature. It used to 
be believed and found on the lattice that {qq) disappeared for T > Tc. However, 
it was found recently, in new accurate calculations for the quenched SU(2) 
case with overlap quarks (not violating CS), that almost a half of the \{qq)\ 
survives at l.bTc >T> Tc- This observation calls for further work both on the 
lattice and in the analytic domain. In addition to the studies with quarks in the 
fundamental representation, discussed above, there appeared several works 
on the phase transition for adjoint quarks, where it was found that < Tsbcs- 
This fact stresses again that we are still far from the detailed understanding of 
the connection between confinement and SBCS. 

(5) An important role in comparison of the QGP properties and results of RHIC 
experiments is played by the transport coefficie nt o f the QGP, and especially 
by the shear viscosity 77 and bulk viscosity C, The existing calculations 
are not very reliable since they exploit nonperturbative dynamics in typically 
Minkowski (not Euclidean) configurations. The FCM, supplied by a careful 
analytic continuation from the Euclidean region, is able to produce both r] 
and C, as saturated by propagators of two- gluon glueballs, coupled by the CM 
confinement. This work is in progress now. 

(6) The FCM approach to the phase transition at nonzero T and /i was given 
above, in section 5.1, based on the assumption that finite density does not af- 
fect strongly the vacuum energy gap, the latter bein g known from conformal 
anomaly. This is supported by lattice calculations of the free energy F{T, fi) 
and looks plausible, since the internal scale of the vacuum energy density is 
connected to excitations with the mass of order of the lowest glueball mass, 
M « 1.5 GeV, whereas the corresponding fi is much lower. However, this ques- 
tion requires further detailed studies since it may change drastically the entire 
picture of the phase transition at low T's due to density. Another point which 
might invalidate the standard picture of the qq pairing phases is the fact that 
vacuum fields, especially measured by the CM correlators, are very strong and 
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can deform or totally destroy this pairing, preferring the qq pairing. 

(7) A very important but still unresolved question is the dimensional transition 
from 4D to 3D with the temperature increase. Studying lattice data for the 
spacial string tension (7h(T) (see, for example, and refs. therein), one can 
see a sharp change of its behaviour at T = 1.5Tc, for SU(2), and at T = 1.2Tc, 
for SU(3), which is not explained by theory. At larger T one has a typical 3D 
behaviour (7h{T) ~ T^. The coefficient was considered, for example, in^^. This 
point of almost "dimensional phase transition" deserves more study. 

(8) There is an interesting development in the theory of the QGP, namely, QGP 
in an external field, which is totally beyond the scope of this review. We refer 
the reader toES and papers cited there for details. 

(9) Finally, let us mention the problem of the phase transition from nuclear matter 
to quark matter. The EoS for the nuclear matter at growing density cannot be 
obtained as an extrapolation of the EoS for the normal nuclear density. As was 
mentioned in , density may strongly distort the usual confinement forces 
between quarks in a nucleon, the latter may loose in its mass and radius, and 
multiquark bags may be formated. This very interesting development may be of 
a practical use and can be connected both to the phenomenon of the so-ca lled 
nuclear scaling and to the hot point of quark stars — see, for example,!^. 

We stop at this point leaving, however, many interesting problems beyond the scope 
of this review. 
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